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INTRODUCTION 
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j 

Growth  in  the  use  of  composite  materials  for  structural  appli-  \ 

cations  has  generated  an  urgent  need  for  the  detection  and  evaluation  of  J 

flaws  Inadvertently  Induced  In  such  structures  by  manufacturing  proce-  ] 

dures.  Such  detection  and  evaluation  procedures  Ccin  be  grouped,  quite 
naturally,  with  other  nondestructive  test  (NDT)  methods  used  in  the 
Inspection  and  evaluation  of  homogeneous  structures.  It  appears  that 
recent  advances  in  the  fields  of  optical  and  acoustical  holography  may 
permit  applications  of  these  technologies  to  the  NDT  field. 

In  recent  years,  filament-wrapped  composite  structures  have  been 
utilized  in  small  solid-propellant  rocket  motors  and  launch  tubes  with 
great  success.  Use  of  composite  materials  In  this  application  has 
merit  due  to  the  high  strenth-to-weight  ratio,  the  low  cost  of  manufac- 
ture, and  the  ease  of  assembly  of  the  rocket  motor  case  and  grain  compo- 
nents. Due  to  the  size  and  nature  of  such  rocket  motors  and  launch 
tubes,  production  would  occur  In  an  automated  fashion  at  very  high  rates. 

Obviously,  such  high  production  rates  will  require  an  automated  system 

of  quality  control.  ' 

It  has  been  found  by  experience  that  flaws  are  likely  to  occur  in 
the  manufacture  of  small  rocket  motor  cases  in  the  form  of  unbonded  sur- 
faces where  bonds  should  be.  These  unbonded  regions  can  occur  for 
several  reasons.  For  example.  If  the  region  is  contaminated  with  a 
foreign  substance  such  as  water,  grease,  etc.,  the  bond  will  not  form. 

Also,  failure  to  apply  adhesive  properly  to  the  entire  region  will 
prevent  the  bond  from  forming.  Such  conditions  are  bound  to  occur 
occasionally  In  a mass-production-type  operation.  The  problem  is  to 
detect  such  flaws,  assess  their  significance,  and  discard  the  faulty 
motor  case  If  conditions  warrant. 

Electro-optical  devices  associated  with  optical  holography  or 
speckle  Interferometry  can  be  used  for  the  detection  of  such  flaws. 

The  optical  holographic  technique  Is  based  upon  the  premise  that  the 
character  of  the  surface  displacements  of  a structure  with  a flaw  are 
different  when  the  structure  Is  loaded  than  when  it  Is  not  loaded. 

Such  differences  In  character  can  be  vividly  detected  interferometrically 
by  comparing  holograms  of  the  loaded  and  unloaded  flawed  structure. 

These  comparisons  produce  fringe  patterns  with  distinct  aberrations 
near  a flaw  [1]. 

At  the  present  time,  the  experimental  techniques  for  detecting 
flaws  have  been  demonstrated  qualitatively.  Assessment  of  flaw  sizes 
and  shapes  have  been  demonstrated. 

The  objective  of  this  work  was  to  determine  the  magnitude  of  the 
surface  displacement  and  surface  strains  of  three  hollow  composite  cylin- 
drical pressure  vessels  in  the  vicinity  of  three  typical  flaws  (or  flaw 
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simulations).  Structural  analysis  to  model  such  structures  is  evalu-  I 

ated  and  suggestions  are  made  for  Improvement.  Section  II  of  this  report 

defines  the  fiber  reinforced  composite  specimens  containing  known  flaws  ■ 

’ including  dimensional  details  and  flaw  sizes. 

The  theory  and  experimental  analysis  of  optical  holography  as  ^ 

applied  to  the  three  specimens  are  Included  in  Section  III.  A technique 
for  real-time  optical  holography  can  be  used  to  observe  how  the  fringe 
pattern  behaves  under  loading  conditions  which  are  time  dependent.  This 
provides  a rapid  way  to  optimize  the  loading  for  a specific  type  of  spec- 
imen. A method  to  determine  the  size  and  shape  of  flaws  from  the  opti- 
cal holographic  fringe  patterns  is  presented.  The  computer  software  for 
this  technique  was  developed  and  a computer  plot  of  "U"  curves  was  used 
to  locate  and  quantify  the  flaws. 

The  equations  needed  to  apply  shearing  speckle  interferometry  to 
thin-walled  cylindrical  launch  tubes  are  derived  in  Section  IV.  This 
method  is  proven  to  be  valid  through  the  use  of  a cantilever  beam  where 
the  displacements  are  easily  measured  or  calculated.  Shearing  speckle 
interferograms  are  presented  for  each  of  three  different  types  of  flaws. 

This  is  a good  whole-field  method  to  detect  and  locate  flaws. 

[ The  theory  for  a single  beam  analysis  using  double-exposure  speckle 

i photography  is  presented  in  Section  V.  Its  validity  is  verified  using 

I a cantilever  beam.  It  can  be  shown  that  this  theory  is  equivalent  to 

I the  point-by-point  data  interpretation  of  Young's  Fringe  Analysis. 

Section  VI  presents  the  theory  for  the  point-by  point  data  inter- 
pretation of  Young's  Fringes.  Equations  to  compute  the  surface  strains 
from  the  fringe  data  for  the  cylindrical  specimen  are  Included.  Sample 
fringe  patterns  and  data  plotted  in  graphic  forms  show  results  of  this 
technique.  It  is  perhaps  the  simplest  method  to  use  to  determine  the 
surface  strains  in  the  region  of  a flaw. 

In  Section  VII,  the  data  reduction  systems  and  software  are  dis- 
cussed in  detail.  System  components  Included  are  the  vldlcon  camera, 
rotating  mirror,  computer  interfaces,  and  the  main  computer  code  descrip- 
tion. The  computer  codes  for  the  vidlcon  camera  and  the  rotating  mirror 
data  reduction  systems  are  Included  in  Appendices  A and  B,  respectively. 

The  system's  hardware  and  operation  including  the  electronic  sub- 
systems are  covered  in  Section  VIII.  Sections  IX  and  X contain  the 

summary,  conclusions,  and  recommendations  for  future  system  improvements.  j 

( 

II.  DEFINITION  OF  THE  PROBLEM  I 

To  aid  in  the  task  of  developing  an  experimental  technique  i 

for  flaw  detection  in  composite  structures,  three  test  structures  were  i 
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designed  with  simulated  flaws  placed  In  specific  locations  and  with 
specific  orientations.  The  simulated  flaws  were  created  by  placing 
thin  strips  of  Teflon  plastic  between  the  two  helical  wraps  of  fila- 
ments at  the  center  of  the  walls  of  the  hollow  cylindrical  structures. 

The  Teflon  strips  were  Intended  to  Inhibit  bonding  of  the  resin  In  the 
composite  between  layers  at  the  location  of  the  strips.  Thus,  the  two 
center  layers  of  the  composite  structure  would  be  free  to  move  relative 
to  each  other.  The  purpose  of  the  three  models  designed  for  study  was 
to  permit  a comparison  to  be  made  of  experimental  and  analytical  pre- 
dictions of  the  displacements  on  the  surface  of  the  cylinders  In  the 
vicinity  of  the  flaws.  All  models  were  manufactured  with  60“  helix 
wraps  and  details  of  each  type  specimen  are  shown  In  Figures  1,  2,  and  3. 


Figure  1.  Composite  cylinder  with  simulated  circumferential  flaw. 
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III. 


OPTICAL  HOLOGRAPHIC  ANALYSIS 


A.  General 

Holographic  interferometry  has  been  used  as  a nondestructive 
testing  method  for  several  years.  The  various  techniques  that  have  been 
developed  essentially  use  a double-exposure  technique  to  produce  inter- 
ference fringes  from  a deformed  body.  If  the  anomaly  in  a body  is  of 
sufficient  size  to  be  critical  to  structural  integrity,  irregularities 
in  the  fringe  pattern  can  be  observed.  These  irregularities  in  the 
fringe  pattern  are  an  indication  of  the  variation  in  the  surface  deforma- 
tion due  to  the  subsurface  void.  Therefore,  if  the  fringes  can  be 
interpreted  to  measure  the  surface  deformation  in  terms  of  displacements, 
some  measure  of  the  structural  integrity  can  be  made  for  the  body. 

Non-destructive  testing  of  thin-wall  cylinders  seems  to  be  a useful 
application  of  holographic  interferometry  because  the  void  must  lie 
close' to  the  surface.  However;  for  filament  wound  composite  cylinders, 
the  interpretation  of  the  fringe  pattern  is  sometimes  complicated  by 
irregularities  produced  by  the  filament  geometry.  In  order  to  determine 
if  holography  could  be  applied  to  the  testing  of  filament -wound  cylin- 
ders, three  test  cylinders  with  three  different  geometries  were  tested. 
Also,  the  displacement  of  the  surface  was  determined  to  measure  the 
surface  variation  in  the  area  of  the  flaw.  Several  photographs  show 
the  different  views  of  each  cylinder  and  illustrate  the  effect  of  flaw 
geometry  and  filament  wrap  angle  on  the  fringe  pattern. 


B.  Application  of  Holographic  Interferometry  to  the  Surface 
Deformation  Measurements  of  Cylindrical  Shells 


1.  Basic  Theory  of  Fringe  Formation 

A common  method  of  forming  double-exposure  holograms 
uses  an  off-axis  reference  beam  which  was  first  developed  by  Leith  and 
Upatnleks  [4].  Other  methods  have  been  developed;  however,  only  the 
off-axis  reference  beam  will  be  considered.  When  a double  exposed 
hologram  is  reconstructed,  interference  fringes  are  formed  on  the  body 
due  to  any  displacement  of  the  body  between  exposures.  If  the  dis- 
placement results  from  deformation  as  shown  in  Figure  4 then  the  fringes 
are  a measure  of  this  deformation.  This  result  is  obtained  due  to  the 
phase  change  of  the  light  due  to  deformation  of  the  body  between  expo- 
sures. The  basic  equation  which  predicts  the  Interference  fringes  is 
the  intensity  of  the  virtual  image.  The  intensity  of  the  virtual  image 
in  a doubly  exposed  hologram  is  given  in  the  following  equation 
[3,  4.  5,  6,  7]: 


^VIRTUAL  = C [1  + cosAe] 


(1) 


where  C is  a constant  and  A9  is  the  phase  change  of  the  light  due  to  the 
deformation  of  the  body, 
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Figure  4.  Geometrical  relationships  for  phase 
change  for  a deformed  body. 


A fringe  will  be  defined  with  From  Equation  (1)  the 

following  condition  results  from  the  definition  of  a fringe 

A©  = (2n-l)jT  n = 1,  2,  3, . (2) 

The  counting  number,  n,  is  called  the  fringe  order.  Equation  (2)  Is  the 
basic  equation  which  relates  the  fringe  information  to  the  deformation 
of  the  body. 

The  phase  change  A © :'.n  Equation  (2)  is  determined  from  the  optical 
path  length  change  of  the  light  between  exposures  of  the  hologram.  In 
the  reference  condition  (undt- formed  state)  the  optical  path  length  is 
denoted  by  Lj^  and  is  the  distance  (Figure  4). 

-4  — ¥ 

= |PS|  + |P0|  (3) 

where  S is  the  light  source 

P is  a point  on  the  surface  of  the  body 
0 is  a point  in  the  hologram  plane. 
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The  distance,  |^|,  between  two  points  (P,S)  is  the  length  of  the  vector 

which  connects  .them.  There£ore>the  distance  |^|  can  be  written  in 
terms  of  the  inner  product  of  the  point  difference.  The  vector  ^ con- 
necting the  points  P and  S is  defined  as  the  point  difference.  Equation 
(3)  can  now  be  written  in  the  following  form: 


= [p^  • pS]  ^ + [p&  • pS] 


1/2 


W 


The  point  P on  the  surface  of  the  body  is  deformed  into  the  point  P' . 
With  the  body  in  the  deformed  configuration  the  optical  path  length  of 
the  light,  denoted  as  L2> is  the  distance 


L2  = Ip'^sI  + \p^o\ 


(5) 


Using  the  inner  product  of  the  point  difference,  this  equation  can 
be  written  as 


p-^s  • P^SJ 


-11/2 


[p^o 


^-11/2 

• p'oj 


(6) 


In  order  to  determine  the  phase  change  of  the  light.  Equations  (5)  and 
(6)  are  used  to  determine  the  change  in  optical  path  length  A,  where 


h -^2 


(7) 


Substitution  of  Equations  (5)  and  (6)  into  Equation  (7)  yields  for 
the  optical  path  length  change 


r_»  _»ii/2  r_ 

A - PS  • PS  + PO  • PO 

} nl/2  r nl/2 

- • pi's]  - [p’>o  • p'>oJ 


1 1/2 


(8) 


The  phase  change  AS  is  related  to  A by  2n/\,  therefore  the  phase  change 
is  written  as 


A©  - ^ I [pS  • pi;]  ^ + [l^  • 

- [p’'s  • - [p^o  • p^o]‘^'|  . 


(9) 
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This  equation,  while  completely  general,  is  not  in  a form  convenient  for 
numerical  use.  Let  the  point  difference  of  the  body  before  and  after 
deformation  be  , then 


p|  = P^S  + P?' 

^ » P^O  + P?'  . (10) 

Therefore,  Equation  (9)  can  be  written  in  the  following  form  using 
equation  (16),  which  yields 


2£  I [pS  + Pt  J ^ + [pt)  • p1)] 


Afl 


|l/2 


- [(P^  • ^)  - 2(pt  •?!')+  (p1>'  • p1")J 

- [(p1)  • pb)  - 2(^  . pt>')  + (pI"  • p1")J 


ni/2 


1/2 


(11) 


In  holographic  interferometry,  displacements  of  points  between 

exposures  are  small;  therefore,  it  will  be  assumed  here  that  jP^'l 

« |P^'|.  Using  this  information  the  last  two  terms  in  Equation  (11) 
can  be  written  in  the  following  form  for  small  deformations  between 
exposures : 


^ ^/2-»  -♦  . p^' 

(PS  • PS)  - 2(PS  . ?P)  PS  • PS — ^ 


PO  • PO 


w 


(12) 


and 


1/2 


^(^  . pb)  -2(pb  . pf ' )J  ==  pb.p1)- 


^ • p?' 


^pb  • pbj 


w 


(13) 


Using  the  results  in  Equations  (12)  and  (13),  the  phase  change  can  now 
be  written  as 


The  tern^  in  brackets  in  Equation  (14)  are  unit  vectors  in  the  direction 
PS  and  PO;  therefore^ the  following  change  in  notation  Is  employed: 


ps 


P|_ 

[pS  • pfe] 


172 


Equation  (2)  can  now  be  written  In  the  following  form; 


(15) 


[ps  + po] 


(2n-l)\ 

2 


(16) 


These  results  are  still  In  a general  form  and  special  cases  will  be 
considered  in  regard  to  the  deformation  of  cylindrical  shells. 


2.  Deformation  of  Cylindrical  Shells 

The  displacement  vector  P^'  of  a point  on  the  surface 
of  the  shell  is  resolved  into  three  components,  a radial  component  a 

circumferential  component  U^,  and  longitudinal  component  U^.  The  experimental 

geometry  and  coordinate  system  of  the  cylinder  are  shown  in  Figure  5. 

The  point  P on  the  surface  can  be  represented  in  either  cartesian  or 
cylindrical  coordinates.  In  the  cartesian  system  let  P?'  have  the 
following  representation: 


P?'  = U e + U "e,  + e. 

1 X2  2 3 

and  the  following  representation  in  the  cylindrical  system 


P?'  = U e + U.  e^  + U "e 
r r $ $ z z 


(17) 


(18) 


The  relation  between  the  unit  base  vectors  are  as  follows: 


-4 

e 

sin 

+ 

cos 

r 

1 

3 

-♦ 

sin 

“ ®i 

cos 

■ 

(19) 
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Figure  5.  Geometry  used  to  measure  the  surface  deformation 
of  cylindrical  shells. 

Let  the  vectors  ^ and  ^ have  the  representations  shown  in  Figure  6. 
Then, 

ps  = ^3®3 

^ + ”“2^2  ^ ”3®3 

where  "*2’  “s  direction  cosines  of  the 

vectors  ^ and  respectively. 

Using  the  notation  for  the  cylindrical  si 
written  in  the  following  form: 

(if  + (U^  sin  $ + cos  + 

+ (i^  + m^)  (U^  cos  ® sin  <I>)  * 
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lell.  Equation  (16)  can  be 


■"2)  U, 


(2n-l)\ 


(21) 


Figure  6.  Geometrical  relationships  for  the  directions 
^ and  vb. 


Equation  (21>  Illustrates  that  the  fringe  Information  In  a single  holo- 
gram contains  Information  for  all  three  displacement  components  that 
cannot  be  determined  uniquely  from  a single  hologram.  In  order  to 
determine  the  components  uniquely  three  holograms  must  be  recorded  simul- 
taneously. Each  hologram  will  have  an  equation  relating  fringe  order 

to  displacements  as  given  In  Equation  (21).  In  each  case,  the  direction  i 

-»  -4  ‘ 

cosines  to  the  unit  vectors  ps  and  po  will  be  different  for  each  holo- 
gram. This  general  procedure  will  require  considerable  data  reduction 

In  order  to  determine  U , U.,  and  U . However,  some  simplifications 

r V z 

can  be  made  In  many  problems;  and  with  careful  selection  of  the  experl-  * 

mental  geometry  the  amount  of  data  reduction  can  be  reduced.  Two  examples 
will  be  considered  to  Illustrate  the  application  of  Equation  (21)  to 
the  measurement  of  surface  deformation  of  cylindrical  shells. 

(a)  Example  1;  Symmetrical  Deformation  of  Cylindrical 
Shells.  The  cylinder  will  be  assumed  to  be  deformed  symmetrically  about 
the  longitudinal  axis.  For  this  case,  the  component  - 0 and  U^, 

are  assumed  to  be  nonzero.  Direction  cosines  for  this  geometry  are 
shown  schematically  In  Figure  6.  The  direction  cosines  for  the  PS 
direction  are 


i 
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Values  for  the  PO  direction  are 
-sin  3 

ni2  * 0 


* cos  3 

This  particular  geometry  assumes  that  the  surface  of  the  cylinder 
is  illuminated  with  collimated  light  parallel  to  the  direction.  The 

direction  of  ^ is  the  same  at  every  point  on  the  surface  of  the  cylin- 
der. Also,  the  direction  cosine  m2  is  assumed  to  be  zero.  This  direction 

cosine  is  the  angle  between  the  observation  direction  and  the  X2  axis. 

If  the  hologram  is  located  at  a long  distance  from  the  cylinder  and  the 
cylinder  axis  X2  is  parallel  to  the  hologram  plate,  then  m2  will  be 

negligible  compared  to  the  other  values.  Therefore,  it  can  be  assumed 
to  be  zero.  The  only  non-zero  displacement  component  to  be  measured  is 
as  shown  schenwitically  in  Figure  7. 


Figure  7.  Synnetrical  deformation  displacement  component  P^' 
for  a circular  cylinder. 
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For  this  example.  Equation  (21)  reduces  to 


-sin  3 + (1  + cos  3)  cos  J ^ ^ 

This  equation  can  be  solved  for  U^,  which  yields 


U * (2n-l)\  o-is 

r 2[(1  + cos  3)  cos  $ - sin  3]  * ' 

Example  1 illustrates  that  for  the  case  of  symmetrical  deformation,  and 
with  propt:r  selection  of  the  experimental  geometry,  U can  be  determined 
from  a single  hologram.  ^ 

(b)  Example  2:  Radial  Deformation  Measurements  of  Cylindri- 
cal Shells  with  Non-Symmetrical  Deformation.  The  second  example  will 
consider  the  general  case  of  the  deformation  of  cylindrical  shells.  Only 
In  this  example  the  radial  deformation  of  cylindrical  is  to  be  measured 
with  non-zero  values  of  U.  and  U . Figure  8 is  a schematic  diagram  of 

9 Z 

the  experimental  configuration  to  measure  the  radial  deformation. 


X,-X2  PLANE 


Figure  8.  Geometry  used  to  measure  the  radial  displacement  component 
for  cylindrical  shells  with  non-symmetrical  deformation. 


For  this  example  the  direction  cosines  are  chosen  to  be  (Figure  9) 


=•  cos  e 


* sin  © 


0 


■ cos  © 


sin  © 


In  this  example  all  light  rays  are  assumed  to  be  collimated  with  Illumi- 
nating and  observation  directions  in  symmetric  directions  with  respect 
to  the  X|  axis. 


•3  4 


' \ 
re  \ 


-Xj  PLANE 


Figure  9.  Geometrical  relationships  for  the  pfe  and  pb 


directions  in  the  - X^  plane. 


Equation  (21)  for  this  example  reduces  to 


2 cos  e (U^  sin  ® cos  9)  - 

If  observations  are  made  In  the  plane,  then  sin  il> 

and  Equation  (24)  reduces  to 

^r  4 cos  6 


(24) 

1,  cos  4 “ 0, 


(25) 


As  In  Example  1,  If  proper  choice  Is  made  with  the  experimental  geometry 
then  the  general  form  of  Equation  (21)  can  be  reduced  to  a form  that  Is 
amenable  to  use  in  numerical  computations.  These  two  examples  are  used 
to  Illustrate  how  the  basic  equation  In  holographic  Interferometry  can 
be  used  to  measure  uniquely  the  radial  component  of  surface  displace- 
ment of  cylindrical  shells. 


A numerical  example  was  chosen  In  order  to  Illustrate  the  application 
of  the  equations  In  Example  2.  This  example  assumes  the  cylinder  to  be 
aluminum  with  the  following  geometrical  and  material  properties: 


1) 

Thickness 

m 

0.020  In. 

2) 

Modulus  of  elasticity 

- 

10.2  X 10^  psl 

3) 

Mean  radius 

■ 

3 in. 

4) 

Poisson's  ratio 

m 

0.332 

5) 

Internal  pressure 

m 

10  pslg 

6) 

\ 

m 

24.9  X 10“^  in. 

7) 

e 

m 

10' 

With  these  properties,  the  fringe  order  can  be  determined  as  a function 
of  the  angular  coordinate.  Figure  10  Is  a graph  of  the  fringe  order 
versus  angular  coordinate.  The  corresponding  fringe  pattern  Is  shown 
In  Section  II. D. 


C.  Surface  Deformation  Measurements  of  Composite  Cylinders 
with  Internal  Flaws 

To  measure  the  surface  deformation  of  composite  Cylin- 
ders with  Internal  flaws,  the  three  test  models  shown  In  Figures  1,  2, 
and  3 were  made  with  simulated  flaws  In  specified  locations.  Each 
cylinder  was  pressurised  and  the  radial  surface  deformation  was 
measured  using  an  experimental  configuration  as  shown  in  Figure  11. 
This  holographic  configuration  will  produce  an  interference  fringe 
pattern  which  Is  a measure  of  the  radial  displacement  of  the  outer 
surface.  The  following  data  were  used  for  each  test: 


9 = 8.5 


X = 632.8  X lo"^  m 

The  measurements  are  made  Independent  of  material  properties.  Only  the 
geometry  and  wavelength  of  light  is  needed  in  order  to  measure  the  sur- 
face deformation. 

The  fringe  pattern  produced  by  a spot  flaw  is  shown  in  Figure  12. 
This  type  of  flaw  is  easily  recognized  from  the  photograph.  The  circular 
fringes  in  the  center  are  approximately  the  same  diameter  as  the  flaw. 


Figure  12.  Fringe  pattern  for  the  cylinder 
with  a spot  flaw. 
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In  the  regions  away  from  the  flaw  the  fringes  are  parallel  to  the  longi- 
tudinal axis.  This  type  of  fringe  pattern  Is  characteristic  of  symmetrical 
deformation.  Therefore,  If  the  cylinder  were  free  of  Internal  voids,  the 
fringes  would  be  a series  of  parallel  lines  away  from  any  end  effects. 

For  this  case  the  radial  deformation  was  calculated  along  a line  In  the 
center  of  the  cylinder.  The  displacements  are  shown  In  graphical  form 
In  Figure  13.  As  can  be  seen  from  this  figure  the  radial  displacement 
changes  by  a large  amount  In  the  region  near  the  flaw. 

Figure  14  shows  the  fringe  pattern  produced  by  the  circumferential 
flaw.  As  can  be  seen  from  the  photograph  this  fringe  pattern  Is  more 
complex  than  the  spot  flaw.  The  Internal  pressure  for  this  cylinder  was 
! 20  pslg.  At  the  higher  pressures  the  deformation  from  the  wrap  angle 

' geometry  becomes  significant.  The  areas  of  the  cylinder  where  the  layers 

overlap  act  as  a stiffener  effect  to  the  shell.  Therefore,  the  fringe 
I pattern  Is  not  one  that  Is  characteristic  of  a symmetrical  deformation  of 

\ an  Isotropic  cylinder.  However,  the  presence  of  the  circumferential 

flaw  can  be  seen  In  the  photograph.  Figure  15  shows  the  radial  deforma- 
tion of  the  cylinder  with  an  Internal  circumferential  flaw.  This  type 
of  flaw  Is  not  as  easily  seen  In  the  photograph.  However,  it  Is  believed 
that  the  spot  flaw  Is  a more  realistic  type  of  flaw  that  will  occur  In  a 
production  situation. 

Figure  16  shows  a composite  cylinder  with  a longitudinal  flaw.  The 
fringe  pattern  for  this  flaw  geometry  Is  significantly  changed  by  the 
presence  of  the  flaw.  This  t3rpe  of  flaw  geometry  produces  a complex 
fringe  pattern  similar  to  the  circumferential  flaw.  However,  this  flaw 
Is  not  expected  to  be  representative  of  a production- type  flaw.  Radial 
deformation  for  this  flaw  geometry  Is  not  as  easy  to  determine  as  In  the  1 

previous  examples.  As  can  be  seen  from  the  photograph,  the  radial 
deformation  Is  essentially  constant  along  the  longitudinal  axis  of  the 
flaw.  Therefore,  the  radial  deformation  was  plotted  as  a function  of 
the  angular  coordinate  relative  to  the  cylinder  centerline.  Figure 
17  shows  the  fringe  order  data  as  It  appears  on  the  photograph.  The 
distance  measured  on  the  photograph  relative  to  the  cylinder  centerline 
Is  then  used  to  determine  the  angular  coordinate.  After  the  angular 
coordinate  Is  known  then  the  radial  deformation  can  be  calculated. 

Figure  18  shows  the  radial  deformation  for  the  longitudinal  flaw  geometry. 

D.  Flaw  Detection  In  Fiber  Glass  Epoxy  Thin-Walled 
Pressure  Vessels 

1.  General 

Weapon  sophistication  points  out  the  need  for  a higher 
strength-to-welght  ratio  and  especially  for  small  structures  ttiat  are 
used  once  and  then  discarded.  This  higher  strength-to-welght  ratio  has 
necessitated  structures  made  from  composite  materials  such  as  glass  and 
epoxy.  However,  along  with  these  more  efficient  structures,  undetermined 
stress  distributions  exist.  The  stress  distribution  solutions  for 
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Figure  16.  Fringe  pattern  for  the  cylinder 
with  a longitudinal  flaw. 
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FLAW  CENTIRUNE 


Internal  pressure  10  pslg. 


homogeneous  structures  are  not  valid.  Until  adequate  analytical  tech- 
niques can  be  developed,  more  testing  is  needed  to  evaluate  structural 
integrity  of  composite -material  structures. 


I 


i 

i 
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For  thin-walled  pressures  vessles  made  from  composite  materials  and 
utilized  as  hand-held  rocket  tubes, safety  considerations  demand  reli- 
ability. Flaws  such  as  broken  fibers  or  delamination  between  the 
glass  and  epoxy  affect  the  structural  Integrity  of  the  structure  and  a 
need  exists  to  detect  these  flaws  that  can  occur  in  manufacturing. 

The  optical  technique  of  holography  is  recommended.  The  technique 
is  nondestructive,  and  can  be  arranged  for  reasonably  fast  evaluation 
on  a real-time  basis.  Further,  the  technique  is  sensitive  enough  to 
detect  flaws  in  thin-walled  structures  without  requiring  pressure  high 
enough  to  degrade  the  structure. 

2.  Holographic  Technique 

In  reference  to  this  problem,  it  was  first  necessary 
to  show  that  flaws  could  be  detected  using  double-exposure  holography. 

A thin-walled  composite  material  cylinder  was  capped  and  pressurized. 
A flaw  was  "manufactured"  in  the  vessel  during  winding  by  Inserting  a 
thin  teflon  disk  (0.5-ln.  In  diameter)  between  the  glass  layers. 

Figure  19  Illustrates  a double-exposed  hologram  of  the  test  vessel.  The 
known  flaw  was  on  the  opposite  side  shown;  thus,  represents  a normal 
vessel  without  flaws.  The  hologram  is  a double-exposed  hologram  with 
one  exposure  taken  at  no  load  and  one  exposure  with  30-psi  load. 

Figure  20  Illustrates  the  presence  of  a flaw  (under  A and  about  the 
center  of  the  vessel). 

The  evidence  of  the  fiber  wrap  angles  between  layers  is  noted. 

When  this  typical  fringe  pattern  is  compared  with  that  of  a typical 
homogeneous  vessel  (Figure  21)  the  difference  is  evident.  While  this 
difference  is  marked  by  careful  restrictions,  some  comparisons  can  be 
made. 
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Observing  Figure  20,  where  the  flaw  is  evident  approximately  half  way 
between  the  end  caps,  the  size  of  the  flaw  (disbond)  is  estimated  at 
0.45-ln.  across.  This  compares  favorably  with  the  0.5-ln.  flaw  manu- 
factured into  the  vessel  (the  diameter  of  the  vessel  is  3.25  in.).  The 
size  of  the  flaw  was  estimated  from  the  photograph  as  the  area  contained 
by  fringes  in  the  vicinity  of  the  flaw.  This  technique  seemed  consistent 
on  other  photographs  of  similar  flaws  and  seems  sufficient  to  estimate 
the  size  of  flaw. 

While  the  technique  of  double-exposure  holography  can  be  used  to 
detect  flaws,  complicating  factors  exist  and  must  be  resolved. 
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One  complication  Is  the  need  to  develop  the  holographic  film,  which 
requires  time  and  thus  Is  not  Ideal  as  an  Inspection  technique  for  mass 
production  of  these  pressure  vessels.  This  complication  can  be  eased 
by  using  a photopolymer  Instead  of  film  which  allows  an  Investigator  to 
observe  the  fringe  pattern  as  It  occurs  (real  time).  One  video  tape 
was  made  of  the  event  to  demons|:rate  capability  and  the  possibility  of 
later  evaluation.  The  preparation  and  use  of  the  photopolymer  material 
Is  described  later  In  Section  V of  this  chapter  with  photographs  of  the 
resulting  fringe  patterns. 

E.  Computerized  Digitization  of  the  Fringe  Patterns 

Another  complication  Is  the  need  to  read  the  fringe  pat- 
terns and  detect  flaws  automatically.  After  considering  several  approaches, 
one  reasonable  way  Is  to  scan  the  fringe  pattern  with  a photometer  at 
small  Increments  along  the  axis  of  the  cylinder  and  note  shifts  from  a 
normal  "s"  curve.  Such  a plot  Is  presented  In  Figure  22.  Line  zero  goes 
through  the  center  of  a flaw  and  Is  located  under  "A"  of  the  figure. 

The  fringe  number  Is  plotted  as  a light  meter  would  detect  a fringe 
scanning  across  the  photograph.  In  a normal  plot  (with  a flaw)  such 
as  along  line  6 the  fringe  spacing  is  close  at  the  beginning,  followed 
by  large  fringe  spacing  and  finally  close  fringe  spacing  as  noted  on 
the  plot.  Evidence  of  a flaw  can  be  detected  on  line  3 where  the  fringe 

plot  deviates  from  the  "s"  curve  and  a "hunqj"  on  the  curve  is  noted. 

This  type  of  effect  is  noted  on  lines  -3,  -2,  -1,  0,  1,  2,  3,  4 and 

established  the  size  of  the  flaw  which  scales  out  to  be  0.59  In.  which 

compares  favorably  since  It  Is  expected  that  the  disbond  Is  slightly 
larger  than  the  known  0.5-in. -diameter  Teflon  in  the  vessel.  Thus,  by 
taking  a photometer  scan  and  digitizing,  the  fringe  spaclngs  (curve 
tangents)  can  be  stored  and  note  made  of  the  deviation  from  the  normal 
"s"  curve  Indicating  a flaw.  In  this  example  the  area  enclosed  by  A. 
represents  the  flaw.  Attention  Is  called  to  a second  trouble  spot 
which  Is  just  below  Aj^  and  labeled  A^.  This  appears  to  be  caused  by 

layers  of  glass  Interacting  at  a bias  angle  and  suggests  a possible 
"flaw"  that  can  be  In  vessels  due  to  the  lay  of  the  windings. 

To  enhance  the  flaw  detection  and  reduce  the  time  required  to  test 
each  cylinder,  computer  software  was  developed  to  scan,  digitize,  and 
analyze  the  holographic  fringe  patterns  in  regions  where  flaws  were 
located  using  a Fringe  Pattern  Analyzer.  The  approximate  size  of  the 
flaws  was  determined.  "U"  curves  were  utilized  in  the  computer  plots 
In  an  effort  to  capitalize  on  the  symmetric  properties  of  the  fringe 
patterns  and  make  flaws  more  easily  detected. 

The  Fringe  Pattern  Analyzer  is  a computerized  system  for  scanning, 
Izlng , processing,  and  displaying  line  Information  recordad  on 
film  ranging  in  size  from  16mm movie  film  to  4 X 5-in.  cut  film  (Flgura  23)  . 

The  scanning  unit  Is  a stepping  spot  scanner  in  which  a spot  Is  stappad 


32 


''VS.- 


over  the  screen  of  a Tektronix  604  display  scope  by  means  ofadlgital- 
to-analog  converter  under  computer  program  control.  The  dlgltal-to- 
analog  converters  have  10  bits  of  resolution  and  drive  the  X and  Y axis 
of  the  display  scope.  The  spot  diameter  of  the  scope  is  approximately 
0.015  in.  which  permits  approximately  300  resolvable  spots  to  be  laid 
down  across  the  5-ln.  dimension  of  the  scope  screen.  The  Image  of  the 
screen  is  then  focused  on  the  film  area  to  be  scanned  and,  conse<]uently, 
the  scanning  resolution  on  the  film  is  approximately  1/300  of  the  long 
dimension  of  the  scope  image  on  the  film.  The  light  which  is  transmitted 
by  the  film  is  collected  by  a large  condensing  lens  system  and  brought 
to  focus  on  a 2-ln.  photomultiplier. 

The  components  of  the  scanning  unit  are  mounted  on  a heavy  2-m 
optical  bench  so  that  they  may  be  accurately  and  reproductively  posi- 
tioned on  the  optical  bench  with  sufficient  accuracy  so  that  no  addi- 
tional focusing  is  necessary.  The  4 X 5-ln.  film  carrier  attachment 
is  loaded  by  placing  the  film  to  be  scanned  between  two  glass  plates. 

Accurate  vertical  and  rotational  positioning  of  the  film  may  be  maintained 
if  the  lower  edge  of  the  film  is  permitted  to  contact  both  of  the  posi- 
tioning pins. 

The  data-collectlon  cycle  is  controlled  by  the  system  software. 

The  data  cycle  commences  upon  Issuing  an  "Intensify"  command.  This* 
command  causes  the  scope  Interface  to  issue  an  "intensify"  pulse  about 
25  nsec  in  duration.  The  light  reaching  the  photocathode  of  the  photo- 
multiplier is  converted  to  an  electrical  signal  and  amplified  Internally 
within  the  tube. 

The  output  of  the  photomultiplier  is  amplified  in  an  amplifier 
mounted  in  the  base  of  the  photomultiplier  carriage.  This  signal  is 
then  sent  to  the  computer  cabinet  in  a coaxial  cable  for  an  additional 
stage  of  amplification. 

After  the  last  amplification,  the  signal  is  fed  to  a sample-and-hold 
circuit  which  is  actuated  by  a "hold"  signal  generated  about  20  ^sec 
after  the  "intensify"  pulse.  The  output  of  the  sample  and  hold  is 
digitized  by  an  8-blt  analog -to -digital  converter  (ADC).  The  conversion 
takes  4 nsec.  When  the  digitizing  is  complete,  the  ADC  sends  a "done" 
signal  to  the  computer  and  the  data  are  transferred  to  the  computer  under  ^ 

program  control.  ' ,> 

j 

The  Fringe  Pattern  Analyzer  is  equipped  with  a Tektronix  4601  hard-  I 

copy  unit  which  prints  out  a graphical  representation  of  the  digitized  \ 

fringe  pattern.  Comparison  of  the  digitized  fringe  pattern  with  the 
original  assures  that  the  correct  fringe  pattern  is  being  analyzed  by 
the  computer. 
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A film  transparency  of  the  fringe  pattern  in  Figure  23  was  digitized 
using  the  Fringe  Pattern  Analyzer.  A computer  graphical  representation 
of  the  fringe  pattern  in  the  region  of  the  flaws  is  shown  in  Figure  24. 
The  scanning  interval  was  0.2  in.  around  the  circumference  of  the  cylin- 
der. The  fringe -density  plots  of  the  "U"  curves  corresponding  to  the 
"s"  curves  plotted  by  hand  are  shown  in  Figure  25.  The  flaws  are  both 


Figure  24.  A computer  graphics  representation  of  the  holographic 
fringe  pattern  from  a composite  cylinder. 
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Figure  25.  Fringe  density  plots  by  computer  of  the  fringe 
patterns  in  Figures  23  and  24. 
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easily  detected  on  the  fringe-density  plots.  With  the  known  scan  inter- 
val and  horizontal  scale  on  the  cylinder  diameter  the  size  of  flaw  can 
be  determined.  Flaw  should  be  circular  and  slightly  larger  than  the 

0.5-ln. -diameter  Teflon  insert  and  is  of  unknown  origin  and  size. 

However,  it  is  nearly  0.5-ln.  and  is  almost  circular.  The  distortion 
in  the  shape  of  the  nearly  circular  flaws  is  due  to  the  curvature  of  the 
cylinder  and  the  unequal  scales  in  the  horizontal  and  vertical  directions. 
It  should  be  noted  that  the  scan  Intervals  for  the  "S"  and  "U"  curves 
are  not  the  same. 

Computerization  of  the  data- reduct ion  process  has  been  demonstrated 
and  this  effort  can  be  extended  to  cover  the  analysis  by  extending  the 
computer  software. 

Other  disbonds  were  manufactured  into  different  pressure  vessels  and 
tested.  For  a disbond  that  occurred  for  several  inches  along  the  length, 
a typical  example  is  shown  in  Figure  26,  The  disbond  is  above  A in  the 
picture.  The  technique  as  described  previously  for  locating  a small 
circular  disbond  will  work  also  in  this  example. 

Another  disbond,  that  of  going  completely  around  the  pressure  vessel 
was  not  easily  detected.  The  point  where  the  Teflon  tape  (under  A and 
about  the  center  of  the  vessel.  Figures  27  and  28)  was  easily  noted 
and  detected  where  the  disbond  was  continuous;  this  was  not  detectable 
with  the  proposed  technique.  However,  it  is  noted  that  such  symmetrical 
flaws  are  not  to  be  expected  and  that  the  beginning  and  end  of  such  flaws 
are  discernable. 


In  an  effort  to  speed  up  the  testing  of  such  vessels  other  load- 
ing techniques  were  examined  to  avoid  the  need  of  special  removable 
plugs  for  pressurizing  the  vessels.  Heating  the  vessel  was  examined  as 
one  alternative,  where  the  disbond  would  cause  different  heat-transfer 
characteristics.  To  date  this  has  not  been  effective.  Vibrating  the 
vessel  is  another  type  of  load  that  was  examined. 


Finally  the  next  phase  of  this  work  is  the  need  to  establish  what 
size  flaw  will  weaken  the  vessels  sufficiently  to  make  them  unsafe. 
This  will  have  to  come  from  testing  to  failure  vessels  that  contain 
various -size  flaws. 

In  summary,  the  technique  of  holography  can  be  used  for  detecting 
flaws.  The  technique  requires  three  views  of  each  vessel  to  include 
120*  in  each  view.  The  results  from  the  double-exposed  hologram  can 
be  recorded  in  real  time  and  stored  on  video  tape.  The  tapes  can  be 
digitized  and  flaws  detected. 
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Figure  28.  Circumferential  disbond  in  a composite  cylinder. 


IV.  SHEARING  SPECKLE  INTERFEROMETRY 


The  theory  and  experimental  results  o£  shearing  speckle  Inter- 
ferometry are  presented.  In  this  section  with  their  potential  for  flaw 
detection  in  composite  cylinders. 

Figure  29  is  a schematic  diagram  of  the  experimental  apparatus  for 
shearing  speckle  interferometry. 


Because  of  the  wedge  in  the  system,  the  following  conditions  will 
be  formed  when  the  lens  plus  the  wedge  images  the  points  P and  P.  on  the 
film  plane.  A shear  in  the  -He  direction  is  assumed; 


a)  The  point  P will  image  at  a point  (x,  y)  and  a point  (X+A.,  y) 
in  the  film  plane. 

b)  The  point  P.  will  image  due  to  the  shear,  at  a point  (x,  y) 
in  the  film  plane. 


Now  the  intensity  of  the  two  images  due  to  double  exposure  of  the 
plate  in  which  the  object  was  deformed  between  exposures  is  considered. 


At  a point  (x,  y)  in  the  film  plane,  the  image  point  of  P and  the 
sheared  i^ge  of  P^  will  have  the  following  light  amplitudes; 


Ep  - A e%^*»  y> 
\ - A e^®p^(*» 


(26) 


The  total  Intensity  of  the  first  exposure  can  be  written  as 


- (Ep  + Ep^)  . (Ep*  + Epj*)  (27) 

where  Ep*  and  Ep^*  denote  complex  conjugates,  or 

- 2A^  + 2A^  cos  0 (28) 

where  0 - 0 (x,  y)  - 0 (x,  y) 

P Pi 

where  now  is  the  intensity  at  the  film  plane  due  to  the  first  exposure. 

Now  the  intensity  for  the  second  exposure  is  considered.  The  light 
amplitudes  for  the  second  exposure  can  be  written  in  the  following  form; 
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To  reveal  surface  deformation,  the  transparency  must  be  spatially 
filtered.  This  is  obtained  by  observing  the  light  amplitude  in  the 
transform  plane.  Figure  30  presents  the  light  amplitude  in  the  transform 
plane  [8], 

2\  « 


G(w)  - e 


■Wj 


g(x)  e“^”*  dx 


(32) 


where  g(x)  ■ a + b I. 


The  following  quantities  are  defined  as 


00 

-/ 


6(w)  ■ I e'"^'”*  dx 


-/ 


C(w)  « t cos  0 e dx 


■/ 


S(w)  ■ f sin  0 e dx 


The  light  amplitude  in  the  transform  plane  can  be  put  in  the  following 
form  / 


G(w)  - e 


-(^)  2 
' ' t(a  + 4bA  ) 5(w)  + 


2A  b C(w)  (1  + cos 


(33) 


+ 2A  b S(w)  sin 


Displacement  information  is  related  to  the  minimum  of  the  light 
intensity  in  the  transform  plane  which  will  be  minimum  when 
(1  + cos  £0)  and  sin  are  both  zero.  Therefore,  when  cos  £0  * 

(2n  - 1)]T,  it  will  define  a fringe.  To  interpret  the  physical  meaning 
of  the  fringes,  first  the  shear  or  deviation  of  the  light  passing  through 
a wedge  is  considered. 

Figure  31  preeeate  the  following: 

» - (eii  - ♦ (e^  - e«) 
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Figure  31,  (/edge  geometry 
The  following  relations  are  noted: 


^312  - I - 9^^ 


A321  - I 


'12 


Jtl32  • a ; 

thus, 

”•(!  - ®tl)  +(f  - »12)  + “ 

or 


a - 0ti  + 0 


and 


i2 


s - «11  + 


t2 


- a 


It  Is  recalled  that 


where  n Is  the  Index  of  refraction  of  the  wedge. 
For  small  angles 


substituting. 


5 - n(0^^  + 0^2)  - “ 


5 ■ (n  - 1)  a 


where 


8 “ the  deviation  angle 
Of  ■ the  wedge  angle 

To  note  the  effect  of  the  wedge  on  the  optical  arrangement,  a 
source  (s)  on  the  lens  axis  that  has  an  image  at  s'  (Thin  Lens  Theory) 
is  considered  (Figure  32) . This  location  was  derived  from 

i + ii  m 2=1 

s s"  (34) 

where  s Is  the  Image  location  for  the  light  source  Intersecting  the 
first  lens  surface  and  assuming  air  as  the  first  medium  and  from 


s"  ^ s' 


for  the  new  wave  surface  leaving  the  second  lens  surface.  Combinina 
Equations  (34)  and  (35) 


7 + 7.  - (n 


(*i  ■ 4)  ■ ' 


where  f Is  the  lens  focal  length.  If  a wedge  Is  positioned  as  Illustrated 
In  Figure  33  the  source  "appears"  to  the  lens  to  be  located  at  s^  instead 

of  s due  to  the  angle  change  (0  - 8)  in  Figure  31.  It  is  also  noted  that 
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*A  “ ®)  " S(0)  for  small  angles 


i + 1,  - J.  + 4 

• • “a  *a 


1 ,.±  + ± 

■ s'  es  s' 

A 

Let  " amount  of  shift,  then  as  can  be  reasoned  from  Figure  33, 

or 

-(f ■']‘®  ' 

Finally, 


- (n  - 1)  os'  (38) 

Next  when  an  object  Is  illuminated  by  a single  beam  as  shown 
schematically  in  Figure  34  due  to  the  wedge,  the  two  points  P and 

are  imaged  at  the  same  point  on  the  film  plane.  The  spatial  separation 
on  the  film  plane  is  denoted  by 

Also  let 

■ s * distance  of  object  (point  P)  from  lens 
■ distance  of  image  plate  from  lens. 

If  the  image  is  not  full  scAle  the  magnification, 

■ ■ D,  A,  , 


can  be  accounted  for  in  the  shift  as 
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(45) 


(n  - 1/2)N  - [(1^  + 1^)  Zat  + ^ 

+ <“8  + “o>  (S  ^ I7  ^ I? 

+ V ^ ^ ^ + I?  ^)] 

It  is  noted  that  a shift  does  not  occur  out  of  the  film  plane. 

An  example  of  a cantilever  beam  with  the  shift  In  the  Z-dlrectlon 
(Figure  36)  Is  considered. 

Fringes  occur  such  that 

(n  - l/2)\  - [(1^  + V ^ + (“a  “o>  S ^“s  + V 1^]  ^ * 

The  direction  cosines  are 

^o-l  * 

1 ■ cos  9 
8 

thus 

(n  - 1/2) \ 

Other  data  Include 

^ - --^2 - 24.91  X 10"®  In. 

cos  9 ■ 0.9421 

sin  0 » 0.3353 

Ab  ■ A ■ (n  - 1)  00 
o'  ' o 

- 0.0899  In. 

u|^  ^ ^ - 0.005  in. 

1 - 6 In. 

V - 0.37 
b - 0.25  in. 


m ■ n ■ o 
o o 


, m ■ o,  Ti  ■ sin  0 
* 8 ’ S 


[(1  + COS  9)^  + sin  0 |j]  As 
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The  displacement  components  are 


U-^  f (1  - Z)(x2  - y2) 

“ H xy^  - (Iz  - J x] 


Let  the  deflection  be  defined  as  6^, 
u (0,  0,  1)  - 5, 


The  displacement  components  now  become 

“ " ^ [2  “ y^>  “ 6 + 2 

- (1*  - 2 *] 

Taking  partial  derivatives 


36,  r 
l3  [2  ( 


1)  (X  - - 


-(1  - s)x] 

Taking  x 0.125  In.  and  y ■ 0 

- -0.200  X lO"^  -34.7  X 10"^  + 416.4  X lo"^  z 

- -51.6  X 10"®  + 8.6  X lo"^  z 


n - ^(7008.8)  + 1210. 


The  theoretical  calculations  are  shown  In  Table  1,  and  a theoretical 
fringe  plot  Is  shown  in  Figure  37.  The  experimental  photograph  Is  Illus- 
trated In  Figure  38.  The  results  were  also  Imaged  using  a TV  system  with 
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a vldlcon  tube;  the  results  are  shown  In  Figure  39.  A comparison  between 
experiment  and  theory  Is  Illustrated  In  Figure  40;  the  agreement  Is  obvious 


TABLE  1.  THEORETICAL  CALCULATION  FOR  CANTILEVER  BEAM 


FIXED  END 


Figure  37.  Theoretical  fringes. 
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CANTILEVER  BEAM  - SHEARING  INTERFEROMETRY 

d -19.59° 

X -6328X 

WEDGE  ANGLE -0.5° 

5 - 0.005  in. 

Figure  38.  Filtered  Image  (all  lines  with  Argon  Laser) 


*11 


Figure  39 


Ftiotograph  of  cantilever  beam  off  TV  monitor 


Comparison  of  theory  and  experiment 


This  technique  was  next  examined  using  a cylinder  (Figure  41). 
The  rectangular  components  of  the  radial  deformation  are 

®x  " \ ^ 

8 " 8 sin  9 

y r 


where 


8 - £2-  A - 

r Et  V 2 ) 


when  the  cylindrical  vessel  Is  assumed  to  be  thln-walled  with  capped 
ends.  In  terms  of  displacement  components. 


“ " - I)  ® 


. 2JL. 


■ 2) 


sin  9 


The  experiment  Is  set  up  as  Illustrated  In  Figure  42.  Unit  vectors  are 
ps  * sin  ^l  + cos  8j 

po  ■ T 


1^  ■ 0 1 ■ sin  P 

o 8 

m_  ■ 1 m <■  cos  6 

o a 

n ■ 0 n ■ 0 

o s 


The  shift  Is  In  the  x direction;  thus  Equation  (45)  Is  expressed  as 

(n  - - ^sln  P ^ + (1  + cos  p)  gj  ^ . < 

The  relation  between  polar  and  cartesian  coordinates  Is  recalled; 


2.2.2 
r ■ X + y 


9 ■ arc  tan  ^ 

X 


and  the  chain  rule  derivatives 
du  du  dr  du 

Sc  Sr  ^ ^ 
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or 


and 


^ ^ sin  d 

^ ^ (cos  6)  - ^ 


I?  ■ ^ (*  • 2 ) “•  ® 


and 


du 


If  " ■ ^ (1  - 2)®^"  ® 

Similarly 


(54) 


(55) 

(56) 


-i) 

sin  9 

and 

■i) 

cos  9 

Equation  (51)  becomes 

L 

1 2/  £x 

Et 

r 

or 

(’■U 

« _EL 

(‘-i 

Et 

r 

(57) 


[sin  p (2x^  + y^)  + (1  + cos  p)  (xy)]  (58) 


+ (1  + cos  p)(x)(r^  - x^) 


1/2 


1 ■ 


(59) 


The  resulting  fringe  pattern  Is  a relatively  complex  function  of  x 
Because  of  this  complexity,  a quantitative  analysis  of  a homogeneous 
cylinder  is  not  straightforward.  In  application  to  a composite  cylinder 
quantitative  analysis  would  be  even  less  feasible.  A qualitative  fringe* 
pattern  for  composites  is  possible  as  noted  in  Figures  43  through  45,  but 
the  complexity  of  assigning  a quantity  to  a given  flaw  makes  this  tech- 
nique less  desirable  than  some  of  the  other  techniques. 
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HORIZONTAL  SHEAR  - WEDGE  ANGLE  0.25° 


0 - 29.67° 

X - 6328  % 

Ap  ° 20  p$i 

Figure  43.  Circumferential  flaw,  filtered  image. 
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HORIZONTAL  SHEAR  - WEDGE  ANGLE  0.25** 
d - 29.67** 

X « 6328  R 
Ap  > 20  psi 

Figure  44.  Longitudinal  flaw,  filtered  image. 


a 

l| 
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HORIZONTAL  SHEAR  - WEDGE  ANGLE  0.25° 


e = 29.67° 


X = 6328  ^ 


Ap  = 20  psi 

Figure  45.  Spot  flaw,  filtered  image . 
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Figure  46.  Arrangement  for 


V.  DOUBLE  EXPOSURE.  SPECKLE  PHOTOGRAPHY 
'SINGLE  BEAM  ANALYSIS  \ 

An  Object  Is  Illuminated  by  a single  beam  as  shown  In  Figure  46. 
The  location  of  the  light  source  Is  shown  schematically  as  S where  P Is 
a point  on  the  surface  of  the  object.  A double-exposure  technique  Is 
considered  vdiere  the  body  Is  Imaged  two  times  on  the  same  film  before  and 
after  deformation.  For  the  first  exposure,  the  complex  light  amplitude 
can  be  expressed  In  the  following  form; 

E (Xj,  X2)  - A (x^,  x^)  e^®(*l'  *2^  (60) 

where  x^^  and  x^  are  the  film  plane  coordinates. 

If  the  object  Is  deformed  and  the  point  P Is  displaced  to  a point 
P , then  the  Image  on  the  film  plane  will  be  shifted  to  new  coordinates 
and  X2»  The  complex  light  amplitude  Is 

E'  (x',  x^)  - A (x^,  xp  (61) 

The  Intensity  of  the  first  exposure  Is 

ij  - i* 


II  ■ A (Xj,  X2) 

The  Intensity  of  the  second  exposure  Is 
Ig  " E'  (I')* 


I,  - x^) 


The  total  Intensity  Is  I where 


IT-I1+I2 


Physically  what  has  been  done  Is  to  Image  the  point  P at  the  location 
Xj,  x^  at  the  film  plane  coordinates.  After  the  body  has  been  deformed, 

the  point  P now  at  P'  Is  Imaged  at  the  film  plane  coordinates  x' , x'. 

-»  12 

If  the  displacement  vector  PP'  is  written  in  the  following  manner 


PP'  • u^  e^ 
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then  the  displacement  vector  in  the  film  plane  can  be  written  as 

[PP'  (x^  x^)]  - M PP'  (x^  X2) 

FILM 

where  M ■ magnification  of  camera,  typically  M < 1 and 

“if  ’ “ “1 
“2f  " ^ *^2 

The  coordinates  at  the  film  plane  for  the  second  exposure  are 
" *1  “if 


(66) 


*2  + “2f 


and 


(PP') 


“if  ®1  “2f  ®2 


(67) 


The  total  intensity  after  the  second  exposure  can  be  written  in 
terms  of  coordinates  in  the  form 


- A (x^,  X2)  + A (x^  + u^^,  X2  + U2P 


(68) 


The  amplitude  transmission  of  the  photographic  plate  is  g(x,  y) . 
Let  the  amplitude  of  transmission ' be  linear  with  exposure  (reasonable 
assumption  for  most  films  that  are  used),  then 


g(Xj^,  X2)  - a + b 
where  a,  b are  constants. 


(69) 


For  convenience  the  x„  coordinate  is  suppressed  in  the  followina 
discussion,  thus 


g(x)  - a + b I A^  (x)  + A^  (x  + 


[' 


“if>] 


Burch  and  Tokarskl  [8]  have  derived  the  following  expression! 
G(w) 


/ a + b ^A^(x)  + A^(x  + Uj^£)J  e^“^'"'^dx 


(70) 


(71) 
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This  expression 
function  of  the 


w 


is  the  Fourier  Transform  of  the  amplitude  transmission 
photographic  plate  where  w is  given  by  (Figure  30); 


e 


2rt  X 

z\ 


(72) 


where  x^  is  the  coordinate  on  the  screen  or  Fourier  Transform  plane  and 

the  limits  of  integration  are  from  - » to  + ».  The  equation  for  G(w)  is 
derived  from  the  Kirchoff  diffraction  formula  and  with  Equation  (70) 

ruwii 

I 2k  I ^ 0 0 

G(w)  - e*-  f a + b A^(w)  + A (x  + Uj^^) 
let 


(-iwx)  , 
e'  ■'  dx 


(73) 


6(w)  ■ / e(-iwx)  dx 

where8(w)  is  called  the  delta  function; 
then. 


+ b e *-  2k  J I ^2^^^  ^ ^2^^  ^ ^ (-Iwx) 

Use  the  shift  theorem  on  the  second  term 

/ [a^  (x  + u^p]  e dx  * 

/ a2  (X)  dx 

Define  F[A]  - / A^  (x)  e dx 

Thus, 


+ b e F[  A]  [ 1 + e ] 

The  intensity  at  the  transform  plane  is 
Ip  ■ G(w)  G(w*) 


(74) 


(75) 


(76) 

(77) 


(78) 


(79) 


or 
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or 


F[A]^  tl  + e 


(-iwUif) 


+ e 


Ip  - 2b^  F[A)^  [1  + COB  (»Ujj)J 
A fringe  will  be  defined  when  I ■ 0;  therefore, 

f 

t 1 + cos  (wUj^p]  ■ 0 

when 

cos  (wUj^P  ■ -1 
or 


(80) 


wUj^^  ■ (2n  - l)jt 
It  is  recalled  that 


(81) 


Figure  47  is  referenced  to  interpret  Equation  (81) . The  information 
at  location  is  selected  by  placing  an  aperture  to  pass  only  the  light 

in  the  immediate  neighbor  of  x^.  This  light  is  then  imaged  on  the  film. 
Thus  x^  being  fixed  the  imaged  fringe  information  is  proportional  to 


where  m^^  is  the  magnification  of  the  new  image  relative  to  the  original 

object.  An  experiment  was  set  up  to  note  the  results  of  the  theory 
reflected  in  Equation  (82) . The  experimental  arrangement  is  shown  in 
Figure  48.  The  coordinates  relative  to  the  cantilever  beam  are  shown  in 
Figure  49. 
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other  data  Include: 


i 

I 


I 

i 

1 


I 


I 


" coordinate  location  on  model  surface 

* - 33  In. 

X "0.98  In. 

8 

\ ■ 20.26  ^lln. 

V - 0.37 

\ 

Scale  factor  ■ 4.55  In.  ■ model  dlmensloiV'photograph  dimension. 
Deflection  on  end  of  cantilever  beam  equaled  0.003  In.  and  0.006  In. 

(In  a direction  parallel  to  film  plane) . Cantilever  beam  dimensions 
were  6-ln.  length,  0.250-ln.  thickness,  and  0.500-ln.  height. 

Experimental  results  are  pictured  In  Figures  50  and  51.  For  compari- 
son with  theory  let 

u^  ■ displacement  In  the  x-dlrectlon 

Uy  ■ displacement  In  the  y-dlrectlon 

u^j^  ■ displacement  In  the  Zj^-dlrectlon 

and  for  the  cantilever  beam  examined 

u , » u ; therefore, 

zl  X 


_ *s  3m6  Tv  ,,  . , 2 2.  1 3 . 1 2 ,1  . 1 

" \*  ^3  U ^ ^ ^ " 6 *1  2 *1  M 2 


(83) 


where  the  displacement  u was  given  by  Equation  (47) . A comparison  Is 

z 

shown  In  Figure  52.  The  agreement  Is  good.  Locating  flaws  In  composite 
structures  by  this  technique  Is  difficult  because  of  the  Indistinct  (rela- 
tive to  holography)  fringes,  but  Is  a full  field  technique. 
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Figure  50.  Cantilever  beam.  Young's  fringes 


Figure  51.  Cantilever  beam.  Young's  fringes 
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VI.  SPECKLE  INTERFEROMETRIC  ANALYSIS  USING  YOUNG'S  FRINGES 


A.  General 

Several  investigators  have  observed  Interference  fringes 
from  laser  produced  speckle  displacements  when  a test  piece  is  loaded 
[9,  10,  11,  12,  13].  Figure  53  illustrates  the  arrangement  for  record- 
ing the  necessary  double-fexposed  film:  one  exposure  with  the  model 
unloaded  and  one  with  the  model  loaded.  Figure  54  Illustrates  how  the 
model  can  be  observed. 

If  each  speckle  point  is  thought  of  as  an  aperture  emitting  light, 
then  several  things  take  place.  First,  the  aperture  (speckle)  produces 
interference  patterns  (Fraunhofer  diffraction)  due  to  its  size;  secondly, 
any  nearby  aperture  will  produce  interference  patterns  with  each  other 
(Young's  fringes).  Finally,  because  the  speckles  (apertures)  are  numer- 
ous and  randomly  oriented  the  patterns  mentioned  previously  are  obscured 
by  each  other.  Now,  if  a film  record  la  made  of  the  speckles  of  a test 
model  and  then  the  model  loaded  and  the  previous  film  record  double 
exposed  (Figure  53),  fringes  can  be  seen  when  coherent  light  is  passed 
through  a small  region  of  the  film  and  the  diffraction  pattern  observed 
(Figure  54).  This  occurs  because  the  film  records  the  aperture  (speckle) 
in  the  undeformed  model  state,  then  the  apertures  (speckles)  are  dis- 
placed in  the  loaded  model  and  again  recorded  on  the  film.  In  a small 
region  of  the  film  the  second  recorded  apertures  are  uniformly  spaced 
relative  to  the  first  set  and  when  a small  coherent  light  beam  is 
passed  through  this  small  region  of  film  a diffraction  pattern  is 
observed  identical  to  those  observed  when  a small  coherent  light  illumi- 
nates a small  region  in  a diffraction  grating.  The  pattern  resulting 
comes  from  interference  of  the  light  as  it  passes  through  a small  aper- 
ture (Fraunhofer  diffraction)  and  interference  from  multiple  apertures 
(diffraction  grating)  uniformly  spaced  in  a small  region. 

The  fringes,  caused  when  light  passes  through  a small  rectangular 
slit  (Fraunhofer  diffraction),  will  initially  be  mathematically  des- 
cribed. Then,  the  pattern  resulting  from  two  apertures  will  be  described 
(Young's  fringes)  followed  by  a mathematical  description  of  the  fringe 
patterns  of  multiple  apertures  (diffraction  grating)  uniformly  spaced. 

B.  Fraunhofer  Diffraction  by  a Rectangular  Aperture 

The  arrangement  shown  in  Figure  55  is  considered.  The 
lens  serves  to  project  onto  a screen  parallel  rays  that  woijld  otherwise 
diverge  from  the  centerline.  A unit  vector  to  the  parallel  light  rays 
along  OH  is  represented  as 


aj^  * il  + mj  + nk  . (84) 


Note  that  OH  Is  on  a ray  passing  through  the  center  of  the  aperture. 
Point  Q Is  a representative  point  In  the  plane  of  the  aperture  and 
located  as 

OQ  = xl  + yj  . (85) 

The  projection  of  OQ  on  to  the  ray  passing  through  the  center  Is 

•-4  — ♦ ^ 

OH  = OQ  * aj^  = xi  + ym  . (86) 


Recall  from  lens  theory  that  the  optical  path  length  QP  and  HP  are 
equal.  From  this,  rays  OP  and  QP  are  out  of  phase  by  CH. 

Let  the  real  part  of 


A Exp  [1  (wt  - ®)]  dydx 


represent  the  optical  disturbance  at  P due  to  secondary  waves  originating 
from  the  element  of  area  dydx  around  the  point  0.  Then,  the  optical 
disturbance  at  P due  to  secondary  waves  originating  at  a small  region 
around  Q Is 


A Exp  [1  (wt  - ^)1  Exp  [ZkI  ^ ^''°)  1 dy  dx 


(87) 


The  Intensity  at  P (being  proportional  to  the  square  of  the  amplitude) 
Is  ° 


= (A  a b)2 

and  the  Intensity  at  P Is 


(88) 


(89) 


To  examine  the  Intensity  on  the  observation  plane  (X,  Y),  note  that  rays 
of  light  passing  through  the  center  of  the  lens  (c)  are  not  bent,  thus 
are  parallel  to  the  original  direction  such  as  a..  The  coordinate  loca- 
tions 5C  and  Y are 


CPj 

ctI 

m « fm 

(90) 


for  small  angles  relative  to  the  center  line  and  where  f Is  the  lens 
focal  length.  Equation  (89)  now  can  be  expressed  as 
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2 


2 


(91) 


Since  the  variables  X and  Y are  separable  In  Equation  (91),  only 


and 


need  be  examined  for  maxima  and  minima  on  the  ^ plane.  Thus,  maxima  and 
minima  occur  when 


sin 

Tf 


and 


0 yielding 


Tf 


= 0,  It,  2n,  — n^^it 


tan  itlU  itXa 

^ ’ Tf 


(92) 


(93) 


From  Equation  (92)  the  initial  value  of  zero  will  produce  a maximum  in 
Equation  (91).  The  other  values  from  Equation  (92),  2it, , nj^’t, 

give  minima  values  of  zero  in  Equation  (91).  The  approximate  values  for 

the  solution  of  Equation  (93)  are  0 (exact),  3it/2 , 5it/2, . All  values 

correspond  to  maxima. 


C.  Fraunhofer  Diffraction  by  Two  Nearby  Apertures 

The  interference  pattern  produced  on  a screen  by  two  equal 
sized  apertures  is  considered.  As  illustrated  in  Figure  56,  consider 
two  apertures  located  along  the  X axis  and  projected  on  the  screen.  Let 
i » sin  ©,  m = 0,  and  n * cos  0;  thus,  the  optical  disturbances  at  P 

from  A,  and  A.,  are 
12 

E|  » A a b \ Exp  [1  (wt  - $,)]  (94) 

\ 


and 


SCREEN 


Figure  56,  Fraunhofer  diffraction  by  two  apertures. 


Exp  [1  (wt  - 


The  phase  difference  at  P due  to  the  different  path  lengths  Aj^P  and 
A^P  is 

2it  |Aj^H|  2n  u sin  © 

* \ 

Then 


^itu  sin  9 


where  I is  the  intensity  at  the  center  of  the  screen  (P  ) «nd  I is 
® _ o 

approximated  as  and  also  equals  sin  6. 

Equation  (96)  contains  the  diffraction  term  as  noted  for  a single 

aperture  and  now  an  Interference  term,  cos^  , due  to  two  aper- 

tures  close  together.  Emphasis  is  called  to  the  cosine  function  where 
zero  Intensity  occurs  if 

- + ^ ~ ^ 

“ 2 sin  B ~ 2x 


,,  . (Is  >1 3 \f 

“ ^2^  sin  © = - f Y 

U S ±( — 2 ) ~ ®i 

X 


= 12  3 


Also,  if  u becomes  less  than  a/2,  then  the  dominating  function  is  the 
sine  function.  The  distance  between  apertures  is  the  Information  that 
will  give  the  displacement;  therefore,  if  the  total  displacement  at  a 
point  is  less  than  one-half  the  aperture  (speckle)  size  then  this  tech- 
nique has  limitations.  A representative  speckle  size  is  300  yin.  and 
150  yin.  should  be  thought  of  as  a minimum  limitation  on  total  displace- 
ment at  a point. 

D.  Interference  Patterns  from  Multiple  Apertures 

Consider  the  final  case  of  multiple  apertures,  the  same 
size  and  equally  spaced  about  the  X-axis  as  illustrated  in  Figure  57. 

The  optical  disturbance  at  P due  to  the  first  aperture  is  given  by 

Equation  (94)  and  for  the  aperture  as 


h ’ Vl  " (wt  - + 2iin^_8ln_©j^  ^ 


The  total  optical  disturbance  from  N diffracted  beams  is 


Ep  * Ej  + E2  + E3  -f  , . . + Ejj 
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and  the  intensity  at  any  point  Is 


I 


P 


(Aab)^ 


- 

- 

, na  sin  6 
sin  7 

A, 

2 

.Njtu  sin  0^ 
sin  ( ^ ) 

irA  sin  6 
\ 

sin 

2 


(100) 


For  the  case  of  U ■ 2a  and  N ■ 10,  an  Illustration  of  Equation  (100) 
Is  shown  In  Figure  58.  It  should  be  noted  that  if  u = a,  no  fringes 
will  occur;  and  represents  the  minimum  limiting  case  for  engineering 
applications.  If  u>a  two  or  more  fringes  will  be  evident  on  the  screen. 
Now  consider  the  application,  if  an  image  is  recorded  of  the  speckle 
pattern  on  a model  made  by  illumination  with  a coherent  light  source, 
these  many  speckles  are  recorded  as  apertures  on  the  film.  When  the 
model  is  loaded,  the  speckles  (apertures)  shift..  By  double  exposing 
the  first  film  record,  a second  set  of  apertures  is  recorded. 

In  a small  region  of  the  film,  the  second  set  of  apertures  is  equally 
spaced  relative  to  the  first.  The  spacing  equals  the  in-plane  displace- 
ment vector.  The  out-of-plane  displacement  effects  only  the  focus  of 
the  image  and  has  almost  no  effect  on  the  aperture  spacing.  When  a 
coherent  light  is  used  to  illuminate,  the  film  record  in  a small  region 
where  the  in-plane  displacements  are  essentially  equal  to  the  projected 
diffraction  pattern  intensity  is  described  by  Equation  (100)  or  as 


itax' 

2 

NnX 

’‘"Tf 

jtaX 

JtXu 

— 

sin  " - 

\f 

(101) 
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Figure  58.  Intensity  distribution  of  multiple  apertures. 

where 


C ■ a proportionality  constant 
a ■ average  speckle  size 

£ ■ lens  focal  length  or  distance  to  the  screen 
\ > wave  length  of  the  light  source 
X « coordinate  distance  on  the  screen 
u In-plane  displacement  . 

Maxima  occur  where 

u - m,^  for  m * 0,  1,  2,  — (102) 

X 

where  Is  the  fringe  order. 
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Instead  of  using  the  lens  focal  length  to  evaluate  the  displacement, 
the  technique  works  without  a lens  and  the  distance  f is  the  distance 
from  the  film  to  the  screen.  This  speckle  technique  was  applied  to  flaw 
detection  in  composite  cylinders.  These  data  are  included  in  Figures 
59  and  60.  The  minimum  spacing  representing  total  displacement  is 
100  pin.  The  maximum  displacements  depend  on  resolution  of  the  fringes 
where  20,000  pin.  is  representative.  Because  sin  0 is  approximated  by 
OT/f) , the  fringe  spacing  is  slightly  nonlinear. 

E.  Fringe  Data  Analysis 

A double^exposed  film  record  was  made  of  the  speckles  on 
the  cylinder.  The  deformation  resulted  from  an  internal  pressure  of 
110  psi.  The  film  was  processed  and  placed  in  a holder  in  front  of  a 
laser  be^un.  A coordinate  system  was  chosen  as  shown  in  Figure  59(b). 

The  fringes  were  projected  on  a screen.  Figure  59(a)  shows  a typical 
fringe  pattern.  The  displacement  was  calculated  in  Equation  (102)  for 
each  location  along  each  line  and  plotted  in  Figure  60. 

The  flaw  location  was  evident  as  shown  in  the  plot  of  Figure  60. 

At  all  lines  within  the  flaw,  the  displacement  was  greater.  Only  two 
lines  are  shown  for  convenience.  In  taking  the  spacing  always  in  the 
axial  direction,  the  normal  strain  in  the  axial  direction  can  be 
calculated  and  used  as  an  excellent  criterion  for  critical  flaws. 

The  experimental  arrangement  is  shown  in  Figure  61.  The  models 
were  loaded  by  pressurizing  the  vessels  to  120  psi.  The  double-exposed 
film  was  analyzed  by  taking  data  at  0.1-ln.  Intervals  (Figure  62).  The 
X,  y axis  was  located  arbitrarily.  Figure  63  shows  the  effect  on  the 
fringes  in  the  vicinity  of  a flaw.  At  each  location  (x,  y) , the  fringe 
spacing  was  measured  in  the  y direction  (axis  of  the  cylinder)  and  dis- 
placement in  the  y direction  calculated  using  Equation  (102).  The 
results  of  the  spot  flaw  illustrated  in  Figure  3 are  shown  in  Figures 
6A  through  69.  The  results  of  the  longitudinal  flaw  Illustrated  in  Fig- 
ure 2 are  shown  in  Figures  70  through  75.  The  results  of  the  circum- 
ferential flaw  Illustrated  in  Figure  1 are  shown  in  Figures  76  through  80^; 

The  results  of  the  spot  flaw  displacements  note  0 to  400  pin.  in 
unflawed  regions.  In  the  vicinity  of  the  flaw  (graphs  in  the  neighbor- 
hood of  X ■ 0.8  through  1.1  in.  and  y ■ 1.1  through  1.6  in.),  the  dis- 
placements change  abruptly  indicating  high  strain.  Figure  63  Illustrated 
how  the  fringes  change.  Figure  69  shows  strain  values  in  the  direction 
of  the  cylinder  axis  (normal  strain  on  an  area  perpendicular  to  the  y 

axis)  that  are  larger  than  2000  (10  ^)  in. /in.  The  flaw  size  and  location 
are  evident.  Similar  results  are  shown  in  Figures  76  and  80  for  the  long- 
itudinal and  circumferential  flaws.  Because  of  the  minimum  stability  (in 
comparison  with  holography  and  correlation  fringe  techniques)  and  because 
of  the  ease  of  quantizing  the  results  in  comparison  with  holography, 
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correlation  fringe,  shearing  speckle  interf eronetry , and  simgle  beam 
speckle  techniques,  this  technique  of  Young's  fringes  is  a most  promising 
method  of  flaw  quantification. 


(a)  TYPICAL  YOUNG'S  FRINGE  PROJECTION 


(b)  GRID  SYSTEM 

gure  59.  Cylinder  with  spot  flaw  and  Young's  fringes 
shown  with  the  grid  system. 


Figure  64.  y-displacement  as  a function  of  y-coordinate  for  the  spot  flaw  cylinder 


Y - COORDINATE  (in.) 

Figure  66.  y-dlsplacement  as  a function  of  y-coordlnate  for  the  spot  flaw  cylinder 
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A -1000<  e < 2000/i-ln./ln. 
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" Figure  69.  Strain  distribution  for  spot  flaw  cylinder. 
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Figure  75.  Strain  distribution  for  longitudinal  flaw  cylinder. 


Y - COORDINATE  (In.) 

Figure  76.  y-displacement  as  a function  of  y-coordinate  for 
circumferential  flaw  (x  - 0,  y - 0) . 
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O -x"O.B 
□ -x-04 
A - * “ 0.7 

V -x-0.8 


Y- COORDINATE  On.) 

Figure  77.  y-d Is placement  as  a function  of  y-coordinate  for 
circumferential  flaw  (x  - 0.5  - x - 0.8). 
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Figure  79.  y-d is placement  as  a function  of  y-coordinate  for 
circumferential  flaw  (*  ■ 1.3  - x ■ 1.4). 
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Figure  80.  Strain  distribution  for  circumferential  flaw  cylinder. 
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VM.  COMPUTER  AIDED  DATA  REDUCTION  SYSTEMS  AND  SOFTWARE 


A . General 

The  computer  aided  data  reduction  system  was  developed 
primarily  for  applications  in  speckle  interferometric  data  analysis.  How- 
ever, other  applications  of  optical  data  reduction  are  amenable  to  this 
type  of  fringe  analysis  such  as  reflecting  grid  and  Moire' [14].  The  basic 
system  consists  of  two  parts  idiich  are  an  optical  recording  of  inter- 
ference data  and  the  computer  hardware  used  in  the  numerical  analysis. 

Optical  recording  of  the  data  can  be  obtained  by  several  methods  depending 
on  the  required  results.  If  only  one-dimensional  displacement  data  using 
speckle  photography  are  desired,  then  a rotating  mirror  system  can  be 
employed.  Two-dimensional  data  analysis  utilizes  a vidlcon  TV  camera 
system  which  is  also  used  in  whole  field  data  analysis  such  as  the  reflect- 
ing grid  and  Moire'.  An  x-y  film  translation  stage  is  common  to  both 
systems  along  with  some  of  the  computer  Interfaces.  This  section  will 
discuss  each  data  recording  system.  Examples  are  selected  which  illus- 
trate the  data  analysis  procedures.  A user's  manual  for  both  systems 
is  also  discussed.  Computer  codes  for  both  systems  are  Included  in  the 
appendices  for  easy  reference.  Applications  of  computer  aided  data 
reduction  to  nondestructive  testing  of  composite  cylinders  are  also 
Included  in  this  section. 

B.  Vidicon  Camera  Data  Recording  System 

A typical  fringe  pattern  used  for  analysis  is  shown  in 
Figure  81.  This  type  of  fringe  pattern  is  produced  when  a laser  beam 
is  passed  through  a photographic  transparency  as  shown  in  Figure  82. 

Interference  fringes  are  a measure  of  the  object  motion  or  surface  defor- 
mation of  the  body.  Theory  for  the  description  of  the  interference 
effects  is  presented  in  Section  VTII.  Because  only  a small  area  is 
Illuminated  by  the  laser,  the  information  in  the  data  analysis  yields  the 
displacement  at  a point.  Therefore,  to  obtain  a complete  map  of  the  sur- 
face, each  point  on  the  film  will  have  to  be  illuminated.  This  is  accom- 
plished by  an  x-y  translation  stage  in  which  the  transparency  has  been  ' 

mounted.  The  stage  has  the  capacity  to  translate  6 in.  in  each  direction  i 

in  0.001-in.  increments  with  a 0.0001-in.  repositioning  accuracy.  Syn- 
chronous stepping  motors  provide  the  control  for  the  translation  direc- 
tions. A photographic  transparency  is  mounted  in  a 4-  X 5-in.  window  as 
shown  in  Figure  83.  The  laser  beam  passes  through  a point  on  the  film 
and  the  film  is  then  translated  relative  to  a stationary  laser  beam.  The 
beam  direction  mirrors  are  attached  to  a post  separate  from  the  trans- 
lation stage.  A photograph  of  the  translation  stage  is  shown  in  Figure  84. 

Location  of  a point  in  the  window  is  controlled  by  the  stepping  motors 
which  are  in  turn  controlled  by  the  computer.  Thus,  the  computer  can  con- 
trol the  translation  distance  with  a minimum  Interval  of  0.001  in.  A 
discussion  of  the  computer  interface  will  follow  in  the  next  section. 
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Light  that  Is  passed  through  the  film  then  forms  a diffraction 
pattern  which  Is  Imaged  by  a vldlcon  camera.  The  Image  on  the  camera 
can  then  be  viewed  on  a TV  monitor  or  stored  In  memory  In  the  computer. 

C.  Rotating  Mirror  Data  Recording  System 

To  Increase  the  speed,  a rotating  mirror  system  Is  utilized 
In  the  data  analysis  of  Interference  fringe  patterns  as  shown  In  Flgure.81. 
A schematic  diagram  of  the  basic  system  Is  shown  In  Figure  85.  This  sys- 
tem differs  from  the  vldlcon  camera  system  In  that  only  a single  direc- 
tion can  be  scanned  by  the  mirror.  However,  the  diffraction  halo  pro- 
jected on  the  screen  (Figure  82)  is  the  same  as  the  halo  projected  on 
the  rotating  mirror  (Figure  85).  As  the  mirror  rotates,  the  fringe 
pattern  Is  swept  across  the  aperture  opening  In  the  photodiode.  The 
aperture  size  Is  sufficiently  small  such  that  only  a small  beam  of 
light  passes  through  the  opening.  Therefore  as  the  mirror  rotates,  the 
fringe  pattern  Is  swept  across  the  face  of  the  photodiode.  The  optical 
performance  of  the  system  was  evaluated  to  determine  some  experimental 
parameters  of  the  rotating  mirror.  The  film  translation  stage  Is  the 
seme  as  discussed  In  the  previous  section. 

The  basic  optical  configuration  of  the  rotating  mirror  analyzer  Is 
shown  In  Figure  86.  A light  ray  from  a divergent  light  source  Impinges 
on  a rotating  mirror  at  an  angle  of  Incidence  )|r  and  Is  reflected  to  a 
photodiode.  The  divergent  light  source  Is  In  practice  the  diffraction 
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ograph  of  the  translation  stage  assembly. 


for  eoaputer  aided  data  reduction 


Figure  86.  Optical  configuration  for  the  rotating  mirror  analyzer. 


halo  emitted  from  a transparency  Illuminated  by  a laser.  The  front  sur- 
faced mirror  rotates  about  an  axis  located  at  a distance  r.  Angle  ^ 

Is  the  rotation  angle  of  the  mirror  and  9 Is  the  angle  between  the  Inci- 
dent light  ray  and  a light  ray  Incident  when  ♦ " 45®.  From  the  geometry 
shown. 


where  the  ip  is  given  as 


"(i) 


(103) 


(104) 


From  the  similar  triangles  AAFR  and  ACER, 

M ^ AF  + CE 
FR  “ AD 


and  it  can  be  shown  that 


AF  -f  CE 
AD 


2A0  cos 


- e)  - 2r  - AC  cos 
AC  sin  ((f>  + v) 


Define  y as 


„ 2A0  cos  ((j)  - e)  - 2r  - AD  cos  (d)  + v) 

AC  sin  (<|>  + v) 

and  letting 

and  substitution  of  M and  N into  the  formula  for  y yields 

(M  -f  /In)  cos  <ti  + (1  + ✓In)  sin  (t  - 2N 
cos  (k  + M sin  6 


(105) 


(106) 


(107) 


(108) 


e » “^  + (^  + tan  ^ y 


(109) 


Figure  87  shows  the  performance  of  a rotating  mirror  for  N ■ 0.25 
and  various  values  of  M.  Results  Indicate  that  for  small  angle  changes 
of  Ip,  angle  6 varies  linearly  as  a function  of  p.  Therefore,  if  the 
divergent  light  field  Is  produced  by  an  Interferogram  and  analyzed  by 
the  photodetector-mlrror  combination,  the  distance  between  maxima  and 
minima  detected  by  the  photodetector  Is  constant. 


The  photodetector  consists  of  a SGD-100  photodiode  whose  output  Is 
amplified  by  an  operational  amplifier  as  shown  in  Figure  88.  The  output 
of  the  circuit  at  v^  Is  then  sent  to  an  analog- to-dlgital  (A/D)  converter 
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Figure  87.  9 versus  ♦ for  N ■ 0.25. 
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Figure  88.  Operational  amplifier. 

for  processing.  Because  the  angle  6 varies  almost  directly  as  a function 
of  (pf  there  Is  very  little  distortion  of  the  electrical  waveform  repre- 
senting the  fringe  pattern. 

D.  Vidlcon  Camera  Computer  Interface 

The  computer  and  Interface  system  Is  shown  schematically 
In  Figure  89.  A PDF  11/40  minicomputer  Is  used  to  control  the  motion  of 
the  translation  table  and  process  the  video  signal.  Operation  of  the 
system  begins  with  a command  from  the  computer  to  move  the  table  to  a 
specified  coordinate  position  In  the  table  window.  With  the  table  In  a 
specified  position,  the  computer  sends  a command  through  the  video  inter- 
face to  store  the  vidlcon  imager  in  the  scan  conversion  memory.  The  scan 
conversion  memory  operates  In  either  a read,  write,  or  erase  mode  depend- 
ing on  the  command  from  the  PDF  11/40.  Upon  receipt  of  an  erase  command, 
the  scan  conversion  memory  switches  to  the  erase  mode  and  the  electron 
beam  is  blanked  and  Is  raster-scanned.  A write  command  Is  then  given  to 
the  scan  conversion  memory  and  the  image  on  the  vidlcon  camera  Is  then 
stored  In  memory.  With  respect  to  the  scan  converter,  read  mode  means 
the  generation  of  an  output  signal  corresponding  to  the  stored  Information. 
The  process  automatically  begins  at  the  end  of  the  write  mode  and  no  read 
command  Is  required.  A command  from  the  computer  then  converts  the  Infor- 
mation stored  In  the  scan  conversion  memory  Into  digital  form.  The  com- 
plete image  Is  not  processed  but  only  selected  lines  of  video  information 
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are  processed  because  of  the  requirements  In  the  data  analysis  which  Is 
discussed  in  the  next  section.  This  process  completes  the  data  acquisi- 
tion at  a point  on  the  film  and  the  computer  then  locates  a new  coordi- 
nate position  and  the  process  Is  repeated. 

Figures  90  and  91  are  photographs  of  the  computer  aided  data  analysis 
system. 


E.  Rotating  Mirror  Analyzer  Computer  Interface 

Figure  92  shows  the  basic  block  schematic  of  the  rotating 
mirror  analyzer.  The  output  of  the  photodetector  (nominal  2 V)  Is  filtered 
with  a low  pass  filter  to  remove  excess  noise.  The  filtered  signal  Is 
monitored  on  an  oscilloscope  and  is  converted  Into  digital  form  by  a 
Digital  AR-11  A/D  converter  whose  output  Is  processed  by  a PDF  11/40  com- 
puter processing  unit.  The  computer  Interfaces  with  a digital  DRll-C 
digital  Input/output  Interface  which  controls  the  scan  converter  Inter- 
face. The  x-y  translation  stage  Is  controlled  by  the  scan  converter 
Interface.  The  motor  assembly  of  the  rotating  mirror  (2-Tln.  diameter) 
rotates  at  39  rpm.  Filtering  by  the  low  pass  filter  Is  accomplished  with 
a 0 to  lOK  ohm  1.0  MFD  reslator-capacltor  circuit.  The  AR-11  digitizes 
data  with  a conversion  time  of  32  pseC. 

In  operation,  the  x-y  translation  stage  Is  advanced  to  a position 
where  data  are  to  be  gathered.  The  rotating  mirror  then  rotates  the 
fringe  pattern  Image  In  front  of  the  photodetector.  The  photodetector 
converts  the  fringe  Intensity  Information  into  electrical  form  which  Is 
then  filtered  to  remove  noise  and  be  digitized  with  the  AR-11  A/D  conver- 
ter. The  digitized  Intensity  Information  Is  then  processed  to  determine 
fringe  spacing  and  localized  displacement  Information.  Because  the  A/D 
converter  and  computer  software  routines  require  a fixed  and  predetermined 
amount  of  time  to  perform  the  digitization  function.  It  Is  not  necessary 
to  use  a line  clock  to  keep  track  of  the  location  of  the  fringe  pattern. 

F.  Vldlcon  Camera  Data  Analysis 

Typical  Information  stored  In  the  scan  conversion  memory 
Is  shown  In  Figure  93  (a).  Data  analysis  consists  of  determining  the 
locus  of  points  of  minimum  Intensity  which  Is  referred  to  as  a fringe. 
However,  this  analysis  is  complicated  by  the  fact  that  the  film  trans- 
parency Is  Illuminated  by  a coherent  light  scarce  which  produces  laser 
speckle.  The  presence  of  this  speckle  Is  to  superimpose  a nolse-llke 
signal  on  the  spatial  Information  which  degrades  the  Image.  In  Fig- 
ure 93  (a),  the  speckle  appears  as  bright  spots  on  the  photograph.  There- 
fore the  effect  of  laser  speckle  is  to  superimpose  bright  spots  on  the 
film  which  complicates  the  location  of  the  fringes. 
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Figure  90.  Video  computer  interface. 


L 

i 


121 


122 


Figure  91.  Complete  system. 


TO-  LOW-PASS 


DIGITIZED  INTENSITY  (V) 


1!^ 


(a) 


SCAN  POINTS  MEASURED  FROM  CENTRAL  SPOT 
(b) 


Figure  93.  (a)  Photograph  of  diffraction  pattern,  (b)  digitized 

scan  of  a horizontal  line. 
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In  the  x-y  read  mode,  the  computer  receives  a voltage  from  the  scan 
conversion  memory  corresponding  to  an  Intensity  at  points  along  a scan 
line.  This  signal  is  stored  in  memory  in  the  PDF 11/40  for  processing 
of  fringe  Information.  A digitized  single  scan  line  of  the  image  In 
Figure  93  (a)  Is  shovm  in  93  (b) . This  signal  Illustrates  the  effect 
of  the  laser  speckle  on  the  Image.  The  Image  should  appear  as  a sine 
wave  without  the  presence  of  speckle;  therefore,  methods  of  reducing 
speckle  noise  are  necessary  to  locate  the  position  of  a fringe. 

A speckle  averaging  technique  was  developed  which  reduces  the 
effect  of  noise  in  locating  points  of  minimum  intensity.  Fringes  which 
appear  as  parallel  lines  on  the  camera  simplify  the  speckle  averaging 
technique.  The  procedure  is  to  digitize  a single  scan  line  and  store 
Into  memory  the  location  and  value  of  intensity  along  this  line.  Another 
scan  line  Is  also  digitized  and  stored  Into  memory.  These  two  scan 
lines,  which  are  separated  by  a small  distance,  are  then  added  and  the 
result  Is  the  average  of  the  two  signals.  The  best  scan  line  separation 
appears  to  be  approximately  10  lines  apart.  The  average  of  these  two 
lines  has  the  effect  of  reducing  the  speckle  noise  superimposed  on  the 
Image.  Figure  94  Is  the  speckle  average  of  two  scan  lines  of  the  Image 
shown  In  Figure  88  (a).  Scan  line  separation  of  this  Image  was  20  lines. 

If  the  scan  lines  of  Figure  93  (b)  and  94  are  compared,  speckle  noise 
Is  reduced  by  averaging  particularly  around  the  points  of  minimum 
intensity. 

Fringe  separation  Is  determined  by  first  locating  the  point  of  maxi- 
mum intensity  from  the  scan  average.  Then  the  first  minimum  point  Is 
located  to  each  side  of  this  maximum  and  the  difference  In  location 
determines  the  fringe  spacing.  The  computer  program  used  for  the  data 
analysis  for  this  system  is  described  In  Appendix  A. 

An  example  of  a cantilever  beam  with  a transverse  end  load  was  chosen 
as  a test  for  the  accuracy  of  fringe  spacing  and  scan  line  separation. 

This  example  was  selected  because  the  theoretical  solution  Is  known 
and  fringe  spacing  varies  from  Infinite  at  the  fixed  end  to  a small  1 

spacing  at  the  free  end.  Therefore,  a wide  range  of  fringe  spacing  >: 

conditions  can  be  examined  from  this  one  example.  ] 

The  beam  geometry  is  discussed  in  Section  IV  and  only  the  results 
will  be  presented  in  this  section.  Figure  95  is  a photograph  of  the 
various  fringe  spacing  conditions  along  the  length  of  the  beam.  Changes 
In  fringe  spacing  along  with  the  presence  of  speckle  can  be  observed  in 
these  photographs.  In  each  photograph  the  horizontal  scan  lines  were 
above  the  dark  spot  In  the  center  of  the  photograph.  Scan  line  separa- 
tion of  10  and  20  were  used  as  a measure  of  speckle  averaging.  Fringe 
spacing  is  determined  by  the  following  [Equation  (102)]  for  " 1 

X - — (110) 

u 
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where 


\ " wavelength  of  light 

f " distance  between  film  and  transform  plane 
u ■ displacement. 


If  the  free  end  deflection  Is  denoted  as  5,  then  the  fringe  spacing  can 
be  written  as 


(111) 


The  ratio  of  displacement  u to  the  free  end  displacement  Is  related  to 
the  fringe  spacing  by  Equation  (111), 


(112) 


This  ratio  was  used  to  determine  the  experimental  results  and  compared 
to  the  analytical  solution.  Figure  96  Is  a graphical  display  of  measured 
values  made  by  an  observer  and  a computer  scan  with  10-llne  separation. 
Both  sets  of  measured  values  are  compared  to  the  analytical  solution. to 
assess  the  accuracy  of  the  measurements.  Compared  to  the  theoretical 
solution,  both  sets  of  measured  values  are  within  2%  of  the  analytical 
solution. 


Results  of  a 10  and  a 20  scan  line  separation  are  shown  In  Figure  97, 
which  shows  that  the  10-llne  scan  yields  the  best  results.  A listing 
of  the  numerical  values  Is  given  In  Table  2. 

G.  Rotating  Mirror  Data  Analysis 


A typical  oscilloscope  trace  of  an  unfiltered  fringe 
pattern  from  the  rotating  mirror  system  is  presented  in  Figure  98,  which 
shows  that  the  presence  of  laser  speckle  produces  a high  frequency  noise 
on  the  signal.  This  signal  Is  filtered  by  a low  pass  filter  to  reduce 
the  speckle  effect  on  the  signal.  Figures  99(a)  and  (b)  Illustrate  the 
comparison  between  a filtered  and  an  unflltered  signal  from  the  fringe 
pattern.  The  filtered  signal  then  served  as  an  analog  Input  to  the 
PDF  11/40  computer.  Fringe  spacing  Is  then  determined  In  this  case  by 
locating  the  maxima  of  the  waveform.  The  maximum  points  are  located  by 
sampling  the  data  and  determining  when  the  slope  between  three  successive 
data  point  samples  changes  from  positive  to  negative.  Appendix  B Is  a 
description  of  the  computational  procedure  used  In  the  data  analysis  of 
the  rotating  mirror. 

To  test  the  accuracy  of  the  rotating  mirror  system,  the  cantilever 
beam  example  was  utilized.  Beam  geometry  Is  discussed  In  Section  IV 
and  only  the  results  will  be  presented  In  this  section.  Fringe  spacing 
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Figure  98.  Oscilloscope  trace  of  the  output  from  the 
photodiode  of  a typical  fringe  pattern  (horizontal 
sweep  5 msec/cm,  vertical  0.5  V/cm) . 


(a)  (b) 

Figure  99.  (a)  Unfiltered  scope  trace  of  the  photodiode  output, 

(b)  Filtered  trace  of  the  same  signal  (horizontal  sweep 
5 msec/cm,  vertical  0.5  V/cm) . 
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ratio  for  this  example  was  determined  by  Equation  (112).  A graphical 
display  of  the  results  is  given  in  Figure  100.  Compared  to  the  theoreti- 
cal solution  the  measured  values  are  within  2%  of  the  theoretical  solution. 


TABLE  2.  RESULTS  OF  SCAN  LINE  SEPARATION  FOR  THE  BEAM  EXAMPLE 


Distance  from 
Fixed  End  Z 
(in.) 

Z/L 

Computer  Spacing 

v* 

20  lines 

10  lines 

20  lines 

10  lines 

3.75 

1 

18 

17 

1.0 

1.0 

3.45 

0.92 

20 

19 

0.90 

0.89 

3.20 

0.85 

23 

22 

0.783 

0.77 

2.95 

0.79 

24 

25 

0.750 

0.68 

2.70 

0.72 

30 

29 

0.60 

0.59 

2.45 

0.65 

32 

35 

0.56 

0.49 

2.20 

0.59 

45 

39 

0.40 

0.44 

1.95 

0.52 

45 

48 

0.40 

0.35 

1.70 

0.45 

64 

61 

0.28 

0.28 

1.45 

0.39 

93 

83 

0.19 

0.20 

1.20 

0.32 

136 

117 

0.13 

0.15 

0.95 

0.25 

- 

190 

- 

0.09 

H.  Main  Computer  Code  Description  ^ 

Two  computer  codes  have  been  developed  to  analyze  inter-  < 

ferograms  using  the  Young's  Fringe  method  of  analysis.  One  of  the  codes 
utilizes  the  vidlcon-scan  converter  while  the  other  controls  the  rotating  , 

mirror  reader  assembly.  These  programs  have  been  written  to  use  the  j 

SYSLIB  library  routines  provided  under  PDF's  RT-11  software. 

Information  for  controlling  the  system  is  passed  from  the  co'^puter 
to  the  scan  converter  interface  via  a DRll-C  digital  input-output  inter- 
face. This  device  resides  at  the  register  address  767772.  Any  infor- 
mation deposited  in  this  location  is  sent  to  the  scan  converter  inter-  \ 

face.  Table  3 indicates  the  available  control  functions  provided.  The 
scan  converter  interface  converts  the  OPCODES  shown  into  commands  which 
it  then  executes  Independently  of  the  computer. 
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TABLE  3.  SCAN  CONVERTER  INTERFACE  CONTROL  FUNCTIONS 


Scan  Converter  Function 

OPCODE 

TV  WRITE 

000001 

TV  ERASE 

000002 

TV  PRIME 

000004 

TV  XY-SEL  READ 

000010 

TV  DC  RESTORE 

000020 

TV  XY  WRITE 

000040 

TV  SELECTIVE  ERASE 

000100 

613  OSCILLOSCOPE  WRITE 

000200 

613  OSCILLOSCOPE  ERASE 

001000 

TV  TRIGGER 

002000 

X-Translatlon  stage  (forward)  direction 
(normal  ■ reverse) 

010000 

X-Translatlon  stage  advance 
(0.001-in.  advance) 

020000 

Y-Translation  stage  (forward)  direction 
(normal  ■ reverse?) 

040000 

Y-Translation  stage  advance 
(0.001- in.  advance) 

100000 

Reset  System 

000000 

Once  a command  Is  sent  to  the  scan  converter  control  Interface,  It 
becomes  necessary  to.  digitize  the  output  of  either  the  scan  converter  or 
the  photodetector  of  the  rotating  mirror  asseird>ly.  This  Is  accom- 
plished by  an  ARll  A/D  converter.  To  start  an  A/D  conversion  for  the 
scan  converter,  the  value  20401  Is  deposited  In  address  770400.  An  A/D 
conversion  for  the  rotating  mirror  system  Is  Initiated  by  depositing  the 
value  21001  In  address  770400.  The  A/D  conversions  are  read  In  address 
770402. 

Six  subroutines  are  used  to  perform  control  functions  for  the  two 
main  programs.  Their  function  and  usage  are  documented  as  follows: 

1)  Subroutine  XADV  (IS,  IR) 

a)  Function  - This  routine  advances  the  X-stage  of  the  XY 
translation  stage. 

b)  Usage  - IS  - number  of  0.001-ln.  Increments  to  be  advanced. 
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(IS  >0)  advanc"  forward 
(IS  < 0)  advance  In  reverse 


IB.  * rate  at  which  the  stage  Is  to  advance. 

10  S IS  i 32000.  The  smaller  values  of  IR  Increase  the  advance  rate. 

2)  Subroutine  YADV  (IS,  IR) 

a)  Function  - This  routine  advances  the  Y-stage  of  the  XY 
translation  stage. 

b)  Usage  - Same  as  for  Subroutine  XADV 

3)  Subroutine  Store  (ISD,  IPD,  IPC) 

a)  Function  - This  routine  primes  the  Hughes  scan  converter 
and  stores  one  video  frame  from  the  vldlcon  camera. 

b)  Usage  - ISD  - duration  of  store  cycle  (4000  ■ optimum  value) 

IPD  - duration  of  prime  cycle  (4000  ■ optimum  value) 
IPC  “ total  number  of  prime  cycles  to  be  performed 

(IPC  > 1) 

4)  Subroutine  Scan  (IX,  lY,  ID) 

a)  Function  - This  routine  digitizes  the  scan  converter 
stored  video  and  stores  the  data  for  further  processing.  The  digitized 
values  have  a range  of  0 to  1024  on  a 0 to  5 V scale. 

b)  Usage  - lY  = locates  the  position  of  a horizontal  scan 
line  (Os  lY  s 1024) 

IX  = locates  the  position  of  a vertical  scan  line 

(0  S IX  S 1024) 

ID  (1,  n)  = horizontal  line  digitized  data  (1  s n S 1024) 

ID  (2,  n)  = vertical  line  digitized  data  (1  S n s 1024) 

5)  Subroutine  PROCES  (ID,  IDD,  IT,  NNX,  NNY) 

a)  Function  - This  routine  searches  the  data  produced  by  sub- 
routine scan  and  locates  the  maxima  and  minima  of  It.  For  example.  If  a 
sinusoidal  signal  is  provided  for  the  filter.  It  will  Icoate  all  the  maxima 
and  minima. 

b)  Usage  - ID  = Data  array  from  subroutine  scan 

IDD  (1,  m)  = locations  of  the  maxima  and  minima  for  the 
horizontal  scan  (Isms  1024) 
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IDD  (2,  in')  = locations  of  the  maxima  and  minima  for  the 
vertical  scan  (1  S m S 1024) 

IT  = maximum  amplitude  of  the  input  signal's  noise 

NNX  s total  number  of  maxima  and  minima  for  the 

horizontal  scan 

NNY  = total  number  of  maxima  and  minima  for  the 

vertical  scan 

6)  Subroutine  RLOAD  (D,  IS) 

a)  Function  - This  routine  converts  the  analog  signal  from 
the  rotating  mirror  scanner  into  digital  form  and  then  calculates  the 
fringe  spacing  from  the  digitized  data  for  the  current  interferogram 
coordinates 

b)  Usage  - D = fringe  spacing 

IS  = trigger  voltage  level  on  a scale  of  0 to  1024 

of  +5  V 

The  ARll  A/D  converter  has  an  input  range  of  0 to  5 V with  0 V 
represented  in  digital  form  as  0 and  5 V as  1024. 

t 

I 

I . Vidlcon  Camera  Code  ! 

The  vidlcon  Young's  Fringe  analyzer  code  is  used  to  con- 
trol the  raster  scanning  of  interferograms  using  the  vidlcon  camera  and 
scan  converter.  The  basic  operations  performed  by  the  program  are  as 
follows: 

1)  Request  scan  parameters  for  the  XY  translation  stage. 

2)  Request  scan  parameters  for  the  scan  converter. 

3)  Request  parameters  for  storing  video  images  and 

filtering  them. 

4)  Positioning  the  XY  table  for  scanning. 

3)  Store  video  Images. 

6)  Digitize  video  images. 

7)  Process  the  digitized  video  for  fringe  spacing  and 
output  the  data. 

Once  the  program  called  VTEST  is  started,  it  requests  the  following 
information: 

1)  X-axis  boundaries,  Y-axis  boundaries  - 415.  This  is  a request 
for  scan  boundary  information.  An  interferogram  is  usually  made  on  a 
4- X 5 -in.  glass  plate  giving  a total  of  4000  X 5000  possible  data  point 
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locations.  Once  the  computer  code  is  started,  the  position  of  the  x-y 
stage  at  that  time  is  assumed  to  be  the  (0,  0)  origin  location.  The 
coordinate  system  for  the  interferogram  is  determined  as  follows.  View- 
ing the  interferogram  along  the  laser  beam  where  it  strikes  the  emulsion, 
the  beam  is  the  z-axis,  the  4-in.  wide  portion  of  the  aperture  is  the 
X axis,  and  the  5-in.  wide  portion  of  the  4-  X 5-in.  aperture  is  the 
Y axis.  The  coordinates  for  a raster  scan  are  input  in  415  format.  For 
example,  to  perform  a full  raster  scan  of  a 4-  X 5-ln.  interferogram, 
with  the  origin  at  the  center  of  the  interferogram,  the  following  data 
are  input:  -2000  2000  -2500  2500. 

2)  X and  Y Increments,  X and  Y axis  stage  rates  - 415.  The  X and  Y 
axis  increments  are  the  number  of  steps  in  the  X and  Y directions  between 
successive  points  where  data  are  collected.  The  X and  Y axis  stage  rates 
control  the  rate  at  vAiich  the  X-Y  stage  advances.  If  in  Step  1,  a data 
scan  was  desired  at  every  other  point  and  the  stages  were  to  move  at  a 
rate  of  10,  the  following  data  would  be  input  to  the  computer;  00002 
00002  0001000010. 

3)  XY-1  and  XY-2  scan  intersects  - 415.  These  parameters  control 
the  position  of  the  scan  lines  for  the  scan  converter.  Two  scans  in  both 
the  horizontal  and  vertical  directions  are  made  to  calculate  fringe 
spacing  and  orientation.  To  scan  at  the  intersection  points  (500,  600) 
and  (510,  610),  the  following  data  are  input;  00500  00600  00510  00610. 

Two  horizontal  scans  at  coordinates  600  and  610  and  two  vertical  scans 

at  coordinates  500  and  510  will  be  made. 

4)  Store  Delay,  Prime  Delay,  No.  Primes,  Filter  level  - 415.  Store 
delay  and  prime  delay,  control  the  erasure,  and  storage  of  a video  frame 
of  data.  Their  input  values  are  recommended  to  be  4000.  The  No.  Primes 
request  is  the  number  of  desired  primes  to  be  performed  between  successive 
video  image  stores.  The  recommended  number  is  3.  Filter  level  is  a 
request  for  how  large  the  amplitude  of  the  noise  is  in  the  signal.  If 

a noise  level  of  0.1  V is  expected,  then  the  computation  of  the  filter 
level  is  as  follows; 

FILTER  LEVEL  - (0.1)  s 20 

Utilizing  the  recommended  input  data, the  data  input  to  the  code  would  be 
04000  04000  00003  00020. 

5)  System  scale  factor  - FlO.4.  This  is  the  scale  factor  for 
converting  the  television  scale  for  measuring  fringe  spacing  into  actual 
physical  units  of  measurement.  To  determine  this  factor  the  computer 
code  may  be  run  initially  with  a scale  factor  of  1.  A test  bar  pattern 
of  known  width  is  then  placed  on  the  diffuser  screen  and  the  computer  . 
then  determines  its  width  in  television  units.  By  knowing  the  fact  that; 
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Localized  displacement 


Scale  Factor 

Television  Spacing  of  Fringes 


the  scale  factor  can  be  determined. 


* 


After  the  scan  parameters  are  input  to  the  computer,  a raster  scan 
will  automatically  occur.  The  switch  console  of  the  PDF  ll/40  coB^>uter 
should  be  at  the  zero  state  at  all  times.  If  Switch  0 is  depressed,  the 
code  will  request  a new  filter  level;  if  Switch  1 is  depressed,  the 
computer  will  rescan  the  same  data  point.  These  switches  are  to  be  used 
if  the  operator  suspects  a problem  In  the  system.  After  each  desired 
point  Is  scanned,  the  coiiq>uter  prints  out  tiie  results.  Depression  of  the 
control  switches  should  be  made  before  the  results  of  the  current  scan 
point  are  completely  printed  out.  A typical  line  of  output  appears 
as  follows: 

X ■ *****  Y ■ *****  DX  ■ *****  DY  ■ *****  XS  “ *****  YS  ■ ***** 
where 

X ■ integer  count  of  current  horizontal  scan  point 
Y ■ integer  count  of  current  vertical  scan  point 
DX  * localized  displacement  In  the  X direction 

DY  ■ localized  displacement  In  the  Y direction 

XS  ■ current  X location  of  the  scan 

YS  ■ current  Y location  of  the  scan  ] 

J.  Rotating  Mirror  Code 

The  second  main  computer  code  Is  for  raster  scanning 
Interferograms  using  a rotating  mirror  asseiif>ly.  The  program  called 
RTEST  requests  the  following  information; 


1)  X>axi8  boundaries,  Y-axls  boundaries  - 415;  (same 
parameters  as  in  VTEST) 

2)  X ''.nd  Y axis  increments,  X and  Y axis  stage  rates  - 
415;  (same  parameters  as  in  VTEST) 

3)  Scan  level,  system  scale  factor  - 15,  FIO.  The  scan 
level  is  computed  from  the  signal  voltage  as  follows;  Measure  the  voltage 
from  the  photodetector  which  clearly  denotes  the  start  of  fringe  infor- 
metion,  then  compute  the  scan  level  from  the  formula; 

Seen  level  ■ Signel  voltage  * 

The  system  scale  factor  is  detensined  tu  in  program  VTEST.  Once  the 
previously  mentioned  three  lines  of  data  are  provided  to  the  computer, 

* rester  scan  will  automatically  occur.  As  In  program  VTEST,  the  console 
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switches  should  all  be  at  their  zero  state.  In  program  RTEST,  depression 
of  Switch  0 causes  a scan  to  be  repeated  while  the  depression  of  Switch  1 
causes  the  computer  to  request  a new  scan  level. 

A typical  line  of  output  data  from  RTEST  appears  as; 

X ■ ****  Y ■ ****  XS  ■ ****  YX  ■ ****  D ■ **** 

where, 

X ■ Integer  count  of  the  current  horizontal  scan  point 

Y ■ Integer  count  of  the  current  vertical  scan  point 

XS  ■ current  X location  of  the  scan 

YS  " current  Y location  of  the  scan 

D ■ localized  displacement  as  determined  by  the  orientation  of  the 
rotating  mirror  analyzer  system. 

The  following  procedure  Is  used  for  running  VTEST: 

1)  Power  up  the  computer,  scan  converter,  laser,  scan  converter 
Interface,  TV  monitor,  vldlcon  camera  and  XY-Table  power  supplies. 

2)  Load  the  RTll  monitor  Into  the  computer. 

3)  Run  program  VTEST. 

4)  Pocus  the  camera  and  position  the  XY-Table  to  Its  origin  point 
for  scanning. 

5)  Place  the  prime,  erase,  TV  write,  XY-Table,  and  TV  XY-SEL  READ 
switches  In  the  comp  position.  Place  all  console  switches  In  their 

0 position. 

6)  Input  the  required  data  to  the  computer. 

The  following  procedure  is  used  for  running  RTEST: 

1)  Power  up  the  computer,  scan  converter.  Interface,  photodiode 
power  supply,  laser  monitor  oscilloscope,  and  X-Y  Table  power  supplies. 

2)  Load  the  RT-11  monitor  Into  the  computer. 

3)  Run  Program  RTEST. 

A>  Position  the  X-Y  Table  to  its  origin  point  for  scanning.  Adjust 
the  signal  filter. 

5)  Place  the  X-Y  Table  switch  to  computer  and  start  the  mirror 
motor.  Place  all  console  switches  In  their  0 state. 

6)  Input  the  required  data  to  the  computer. 
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K.  Applications  to  Nondestructive  Testing  of  Composite 
Cylinders 

A computer  scan  was  made  of  the  speckle  photograph  of 
composite  cylinders  with  simulated  flaws  shown  in  Figures  2 and  3. 

Speckle  photographs  of  the  cylinders  were  the  same  photographs  used  in 
the  data  analysis  presented  in  Section  VI.  Computer  scans  were  made 
along  several  sections  of  each  cylinder  to  determine  if  the  data  analysis 
was  sensitive  to  the  simulated  flaws. 

Figures  101  and  102  are  graphical  representations  of  scan  lines  from 
the  rotating  mirror  system.  Scan  line  directions  for  this  application 
was  along  the  X direction  as  shown  in  Figure  62.  Figure  101  presents 
the  results  of  the  fringe  spacing  for  the  simulated  spot  flaw.  Figure  74 
presents  the  results  of  fringe  spacing  for  the  longitudinal  flaw.  These 
two  figures  show  that  the  computer  scan  is  capable  of  detecting  the 
flaws  to  the  same  accuracy  as  an  observer  recording  the  data. 


VIII.  COMPUTER  AIDED  DATA  REDUCTION  SYSTEMS  HARDWARE  AND 

OPERATION 

A.  General 

The  purpose  of  this  section  is  to  provide  technical 
information  about  the  video-computer  interface.  Model  4000,  and  to  aid 
technical  personnel  in  the  operation  and  understanding  of  its  electronic 
circuits . 

This  section  contains  a general  description  of  the  Model  4000  and 
the  speckle  analyzer  system.  A detailed  description  of  the  operation 
of  each  major  circuit  within  the  Model  4000  Interface  is  also  provided. 

All  details  needed  to  fabricate  or  purchase  the  components  of  this 
system  have  been  completely  documented  but  only  the  block  diagrams  and 
control  panel  layouts  are  Included  in  this  report.  System  schematics, 
test  points,  cable  systems,  interfacings,  pin  connections,  and  catalog 
data  on  other  components  are  included  in  Reference  15.  Complete  fabri- 
cation drawings  were  prepared  and  can  be  made  available  if  components 
of  this  system  are  needed. 

B.  System  Description 

The  video-computer  interface  panel.  Model  4000,  was  built 
for  the  Ground  Equipment  and  Missile  Structures  (GE&MS),  Directorate  of 
MIRADCOM  by  Sperry  Support  Services  (SSS) . The  system  consists  of  the 
following  items: 


Flgtire  102.  Rotating  mirror  scan  results  of  the  longitudinal  flaw, 
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Supplier* 

1)  Video-computer  interface  panel,  A 

Model  4000. 

2)  X-Y  translation  table.  A 

3)  Television  camera,  optical  acessories,  A,B 

STV  1000. 

4)  Minicomputer,  POP  11/40.  B 

5)  Laser,  Model  124A.  C 

6)  Scan  converter.  Model  639.  D 

7)  Video  storage  unit.  Model  613.  E 

Supplier's  Code:  (A)  SSS 

(B)  MIRADCOM,  GE&MS  Directorate 

(C)  Spectra  Physics,  Inc. 

(D)  Hughes 

(E)  Tektronix,  Inc. 

The  video-computer  Interface  is  shown,  as  it  fits  into  the  Electro- 
mechanical Single-Beam  Speckle  Interferometric  Analyzer  System,  in 
Figure  103.  All  the  separate  pieces  of  equipment,  purchased  by  different 
companies,  are  tied  together  to  the  PDP  11/40  minicomputer  with  the 
Model  4000  interface. 


The  laser  passes  a beam  of  light  through  a hologram  specimen,  and 
is  reflected  onto  the  diffusion  screen.  The  vidicon  television  camera 
makes  a picture  of  the  image  on  the  screen.  The  scan  converter  then 
stores  the  picture.  A horizontal  and  vertical  line  is  scanned  and 
digitized  by  the  minicomputer.  The  results  are  processed  and  read  out 
on  the  terminal.  The  stored  video  can  be  seen  on  the  monitor. 

C.  Electronic  Subsystem  Functions 

The  purpose  of  the  video  computer  interface  is  to  provide 
the  required  interfacing  between  a PDP  11/40  minicomputer  and  other 
equipment  of  the  system  shown  in  Figure  103.  The  interface  also  provides 
for  the  manual  operation  of  all  the  additional  equipment  from  a single 
source . 

Included  within  the  video  computer  interface  panel  is  a video 
amplifier,  a video  filter,  a video  comparator,  a board  of  electronics 
for  manual  and  computer  commands,  and  a board  for  motor  drive  commands. 
The  block  diagram  of  the  subsystems  is  shown  in  Figure  104, 
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Figure  104.  Video  computer  interface  block  diagram. 


Also  provided  is  an  X-Y  translation  table  and  power  supply.  The 
X-Y  table  is  capable  of  moving  a specimen  in  two  dimensions,  covering 
2 

100  in.  , within  0.001  in.  in  either  direction.  The  power  supply  is 
located  on  a separate  panel.  It  provides  the  current  to  the  drive 
motors . 


D.  Circuit  Description  of  the  Video  Computer  Interface 

The  electrical  block  diagram  for  the  Sperry  Model  4000 
video  computer  interface  is  shown  in  Reference  15 . Also  shown  are  its 
interconnections  to  the  X-Y  table,  TV  camera,  PDF -11/40,  scan  converter, 
monitor,  and  storage  monitor. 

The  motor  drive  logic  card  provides  the  axis  control  and  axis  direc- 
tion logic.  The  control  commands  may  originate  from  either  the  mini- 
computer or  from  the  front  panel.  The  outputs  go  to  the  X-Y  Translator 
Driver. 

The  command  logic  card  gates  the  commands  from  the  computer  or  from 
the  front  panel  switches  to  control  the  other  subsystems.  The  subsystems 
controlled  by  this  card  are:  Scan  Converter,  Storage  Monitor,  and  TV 
camera.  The  light  emitting  diodes  (LED)  drivers  are  located  on  the  card 
also. 


The  video  comparator  card  is  a video  peak  detector.  The  output 
is  a coded  digital  signal  resulting  from  the  comparison  of  real  time 
video  to  a reference  voltage  level. 

The  two  video  amplifier  cards  provide  for  the  polarity  and  amplifi- 
cation of  the  input  video.  The  amplifiers  provide  a 75-ohm  terminated 
input  and  are  capable  of  driving  a 75-ohm  output.  The  amplifiers  provide 
a variable  gain  between  0.3X  to  5X. 

The  video  filter  card  provides  a variable  bandpass  filter.  The  center 
frequency  is  variable  from  500  kHz  to  8 MHz;  the  bandwidth  is  variable 
from  200  kHz  to  8 MHz  (2  pole  12  dB/octaT-e) . Also  available  on  this  card 
is  the  sync  stripper  and  sync  clamp  to  remove  the  sync  from  the  output 
video. 

The  translator  module  and  servo  motors  are  purchased  items  from 
Superior  Electric  Company.  Instructions  for  operation  and  maintenance 
may  be  found  in  Reference  15. 

E.  Theory  of  Operation 

1.  Motor  Drive  Logic 

The  motor  drive  card  provides  the  gated  clocks  to 
the  translator.  These  clocks  are  called  X-Motor  Drive  (XMDOO,  Forward; 


XMDOl,  Reverse)  and  Y-Motor  Drive  (YMDOO,  Forward;  YMDOl,  Reverse).  They 
are  generated  from  XMLOO  and  YMLOO,  X-Motor  Logic  and  Y-Motor  Logic. 

The  direction  of  the  motors  is  determined  by  the  controls  on  the  front 
panel . These  commands  are  X-Direction  (XDSOO)  and  Y-Direction  (YDSOO) . 
The  X manual  drive  commands  XASOO  and  XAROON,  set  the  S-R  flip-flop  Al. 
The  Y manual  drive  commands,  YASOO  and  YAROON,  set  a similar  Hip-flop 
Al.  The  outputs  of  the  flip-flops  are  further  gated  with  the  Computer/Manual 
signal  CMSOO  to  prevent  the  manual  operation  during  the  computer  control 
period.  The  computer  commands  are  similar  to  the  manual  commands,  but 
are  useable  only  when  the  COMPUTER/^NUAL  M)DE  switch  is  in  the  COMPUTER 
position.  These  commands  and  their  functions  are  as  follows;  XSAOO 
and  YSAOO  are  the  X-Y  stage  advance;  XSDOO  and  YSDOO  are  the  stage 
direction. 

2.  Command  Logic 

The  command  logic  card  ties  the  Hughes  Scan  Converter 
and  the  Sperry  STV-1000  camera  to  the  minicomputer.  The  manual  operation 
of  the  equipment  is  possible  through  this  card.  The  prime,  erase,  write 
outputs  to  the  scan  converter  are  generated  from  the  computer  and  from 
the  front  panel  switches.  The  computer  commands  are,  respectively,  CTPOO, 
CTEOO,  and  CTWOO.  The  manual  commands  are  TPCOO,  TECOO,  TWCOO.  These 
inputs  are  "OR' ED"  together  by  the  diodes  CRl  through  CrH  The  resulting 
outputs  then  are  HTPOO,  HTEOO,  and  HTWOO,  respectively.  The  dc  restore 
may  also  be  generated  by  the  computer  or  manual  switches.  These  commands 
trigger  the  one  shot  A3  to  provide  a 40-n8ec  pulse.  The  output  to  the 
scan  converter  is  HRCOO.  There  is  a provision  for  manually  erasing  and 
writing.  The  controls  are  located  on  the  front  panel.  The  inputs  to 
this  card  from  these  switches  are  called  ERSOO  and  WRSOO.  Extreme  care 
should  be  taken  with  the  MANUAL  WRITE  switch.  Should  this  switch  be 
activated  for  more  than  an  instant,  permanent  damage  may  occur  to  the 
scan  converter.  The  LEDs  are  driven  by  A2,  A4,  and  A5.  The  LAMP  TEST 
switch  on  the  front  panel  lights  all  the  LEDs  on  the  video  computer 
Interface  panel. 

3.  Video  Comparator 

Video  peak  detection  is  accomplished  with  the  use  of 
this  card.  The  live  video  is  brought  into  Comparator  Al.  This  signal 
is  compared  with  the  threshold  voltage  set  by  Potentiometer  R3.  The 
analog  ground  is  kept  isolated  from  the  digital  ground  on  this  card  by 
the  use  of  Optoisolator  A2.  The  output  signal  of  A2  has  a slow  rise-time 
and  a slow  fall-time  due  to  its  Inherent  characteristics.  The  signal 
is  fed  into  a Schmlt-trigger  Inverter  to  provide  a better  digital  repre- 
sentation of  the  video  signal. 

4.  Video  Amplifier 

The  video  amplifier  card  provides  a voltage  gain  of 
2.  This  gain  is  provided  by  Video  Amplifier  Al.  The  output  polarity 
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may  be  changed  by  the  POLARITY  switch  on  the  front  panel.  This  switch 
applies  the  video  signal  to  the  Inverting  or  noninverting  inputs.  A2 
is  a high  speed  buffer  amplifier.  It  allows  the  video  amplifier  to  drive 
75-ohm  cables. 

5.  Video  Amplifier  2 

This  card  is  identical  to  the  Video  Amplifier  1 with 
a few  exceptions.  The  first  of  these  is  that  the  differential  Inputs 
are  not  available  to  Al.  Secondly,  the  overall  gain  of  five  required 
by  the  system  is  achieved  by  this  card. 

6.  Video  Filter 

The  video  filter  is  made  up  of  a high  pass  filter 
and  a low  pass  filter  in  series.  The  low  pass  filter  is  composed  of  Video 
Amplifier  Al  and  the  resistor-capacitor  (RC)  feedback  network  as  selected 
by  the  UPPER  CUT  OFF  FREQUENCY  controls.  The  high  pass  filter  is  composed 
of  Video  Amplifier  A2  and  the  RC  feedback  network  as  selected  by  the 
LOWER  CUT  OFF  FREQUENCY  controls.  The  sync  is  removed  by  Comparator 
A3  and  then  it  is  used  to  gate  out  the  sync  in  the  video  before  sending 
it  to  the  video  filter. 

F.  Electronic  Subsystem  Controls 

Before  activating  the  video  computer  Interface,  the 
following  initial  conditions  must  exist: 

1)  All  electronic  cards  are  inserted  in  their  proper 

positions . 

2)  All  cables  are  run  correctly  with  their  connectors 

securely  mated. 

The  video  computer  Interface  panels  are  shown  in  Figures  105  and  106. 

1.  Operating  Control  Functions 

The  POWER  switch  on  the  Model  4000  provides  the 
primary  power  to  the  ±15  V and  +5  V power  supplies.  When  this  switch 
is  in  the  ON  position,  the  115  V light  as  well  as  the  +15  V,  -15  V,  and 
+5  V LEDs  should  turn  on.  < 

The  MANUAL  TABLE  CONTROL  portion  contains  the  manual  controls  to 
operate  the  X-Y  Table.  These  controls  are  useable  only  when  the  TABLE 
CONTROL  switch  is  in  the  MANUAL  position.  The  available  controls  are 
X and  Y directions  (FWD  and  REV),  SPEED,  and  ADVANCE.  The  INDEX  LEDs 
turn  on  when  the  X-Y  Table  has  returned  to  its  home  position. 
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The  SCAN  CONVERTER  section  contains  the  switches  to  activate  the 
scan  converter  from  the  computer  or  manually.  The  operations  available 
are  WRITE  and  ERASE.  Extreme  care  should  be  taken  to  prevent  the  manual 
writing  of  the  scan  converter  for  longer  than  a few  milliseconds  in  one 
location.  Permanent  damage  may  result  if  this  switch  is  left  in  the 
MANUAL  position  for  only  a few  seconds.  The  PRIME  command  provides  a 
prime  signal  and  an  erase  signal  to  the  scan  converter.  The  X-Y  READ 
may  also  be  controlled  by  the  computer  and  by  the  manual  switch.  The 
direct  current  RESTORE  pulse  is  generated  by  the  push  button  on  the 
panel. 


The  613  SOURCE  portion  of  the  video  computer  interface  controls  the 
Tektronix  613  storage  monitor.  The  613  can  be  controlled  in  three  ways. 
These  are*,  the  minicomputer,  the  manual  switches,  or  by  the  video 
representation  of  the  video  comparator.  The  two  functions  controlled  by 
these  controls  are  WRITE  and  ERASE. 

The  CAMERA  portion  allows  the  STV-1000  TV  camera  to  be  operated 
by  the  computer.  TRIGGER,  EXTERNAL  BLANKING,  and  FLASH  may  be  done 
externally  from  the  camera. 

The  video  filter  section  contains  the  controls  to  adjust  the  band- 
width of  the  video  signal.  The  two  switches  to  select  the  cut  off  fre- 
quencies and  the  two  variable  controls  allow  for  a bandwidth  from  0.5 
to  8 MHz  or  any  bandwidth  in  between.  These  controls  are  labeled  UPPER 
and  LOWER  CUTOFF  FREQUENCY,  The  VIDEO  INPUT  switch  provides  for  the 
selection  of  two  different  video  signals  A and  B.  The  PROCESSOR  GAIN 
determines  the  system  gain  CO. 5 to  5X) . The  MONITOR  switch  selects  the 
type  of  output.  The  INPUT  LOOP  mode  bypasses  all  the  electronics  and 
sends  the  live  video  straight  to  the  monitor.  The  PROCESS  OUT  mode  allows 
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the  video  filter  (and  also  the  video  gain)  to  be  switched  In.  Positive 
( and  negative  outputs  from  the  scan  converter  can  be  selected  by  the 

[ SCAN  CONVERTER  VIDEO  POSITIVE/NEGATIVE  switch. 

2.  Alignment  Procedures 

All  electronic  subsystem  circuits  are  designed  to 
operate  without  the  need  for  any  special  adjustments  beyond  those  pro- 
vided by  the  front  panel  controls.  However,  the  threshold  to  the  video 
comparator  card  (A4)  may  be  adjusted  to  provide  a more  sensitive  digital 
representation  of  the  video. 

To  adjust  the  threshold,  a small  screwdriver  should  be  used  to  turn 
the  potentiometer  on  Card  A4.  Turning  the  screwdriver  counterclockwise 
causes  less  sensitivity  to  light.  At  the  upper  limit,  the  output  will 
correspond  to  all  black;  the  lower  limit  will  correspond  to  all  white. 
Proper  adjustment  will  cause  a pulse  at  Test  Point  1 that  corresponds 
to  a peak  in  the  video.  The  voltage  should  be  approximately  0.5  V at 
Pin  2 of  Al. 

IX.  SUMMARY  AND  CONCLUSIONS 

Approximately  thirty  fiber  reinforced  composite  cylinders  were 
fabricated  in  the  in-house  composite  structures  facility.  Three  types 
of  programmed  flaws  were  fabricated  into  the  composite  structure.  The 
flaws  consisted  of  delaminations  between  the  fibers  and  were  produced 
by  inserting  some  Teflon  tape  in  a 0.04-in. -thick  cylinder  wall.  This 
type  of  flaw  corresponds  to  those  resulting  from  dirt,  grease,  voids,  or 
flaws  being  fabricated  Into  the  structure.  Spot,  circumferential,  and 
longitudinal  flaws  were  used. 

The  basic  theory  for  the  analysis  of  the  optical  holographic  fringe 
patterns  was  adapted  to  the  missile  launch,  tube  or  motor  case  geome cries 
and  boundary  conditions  for  the  specimen  Involved.  Equations  to  reduce 
the  fringe  patterns  to  surface  displacements  were  derived.  Test  fixtures 
and  apparatus  were  desl^nf'^  and  fabricated.  Optical  holographic  fringe 
data  were  taken  from  specimens  of  each  type  with  pressure  loadings  using 
both  the  double  exposure  and  real-time  methods.  Fringe  patterns  were 
so  sensitive  that  at  certain  loadings  even  the  wrap  angles  of  the  fibers 
and  the  crossover  regions  were  easily  detected.  Because  locations  of 
these  points  are  well  known,  interference  with  the  detection  of  flaws 
is  not  considered  a major  problem.  The  real-time  approach  is  a rapid 
way  of  optimizing  the  loading  parameters  to  best  resolve  the  flaws. 

Each  type  of  flaw  was  detectable  using  double-exposure  holography  and 
real-time;  however,  both  methods  require  careful  digitization  and  analysis 
to  quantify  the  flaws  and  determine  if  they  are  critical. 
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Shearing  speckle  and  single  beam  speckle  Interferometry  have  both 
been  demonstrated  to  be  able  to  detect  the  presence  of  flaws  In  fiber 
wound  composite  cylinders.  Although  each  can  detect  flaws  the  two 
experimental  techniques  are  best  adapted  to  separate  functions  because 
they  complement  each  other.  Therefore,  In  an  overall  nondestructive 
testing  program  each  technique  can  be  utilized  to  perform  specific  func- 
tions . 

The  basic  theory  and  applications  of  speckle  shearing  Interferometry 
are  presented  In  detail.  This  technique  has  the  ability  to  detect  flaws 
from  a whole  field  approach  with  good  fringe  visibility.  Furthermore, 
with  a combination  of  mirrors  the  complete  cylinder  can  be  observed  and 
photographed  on  a single  photographic  record.  Thus,  a rapid  visual 
observation  can  be  made  for  Initial  Inspection  of  Internal  voids.  The 
analysis  of  fringe  formation  to  the  deformation  of  circular  cylinders  was 
also  presented.  The  resulting  fringe  pattern  Is  a relatively  complex 
function  of  the  cylinder  geometry  and  assigning  fringe  orders  Is  not  as 
convenient  for  computer  aided  data  reduction.  However  qualltlve  analysis 
Is  very  feasible  using  this  technique. 

A complete  discussion  of  the  single  beam  speckle  photography  Is 
Included.  Double-exposure  speckle  photography  is  a whole  field  observa- 
tion and  the  Young's  Fringe  technique  is  the  point  by  point  data  inter- 
pretation. Results  of  the  speckle  photographic  technique  Illustrate  that 
Internal  voids  can  be  detected  and  quantified  In  terms  of  the  surface 
strains  of  the  cylinder.  Data  analysis  Is  the  easiest  to  apply  computer 
aided  data  reduction.  Equations  are  easy  to  Interpret  and  only  fringe 
spacing  Is  required.  Fringe  orders  do  not  have  to  be  assigned  thus 
allowing  for  straightforward  data  Interpretation.  This  technique  repre- 
sents a method  to  quantify  the  flaw  and  assess  the  structural  Integrity 
of  the  cylinders  with  the  presence  of  Internal  voids.  Data  were  presented 
to  validate  each  technique. 

A complete  nondestructive  testing  system  to  detect,  locate,  size, 
and  quantify  flaws  In  fiber  reinforced  composite  missile  launch  tubes 
and  motor  cases  was  designed  Including  both  the  hardware  and  software. 

The  system  was  assembled.  Interfaced,  and  checked  out.  Most  of  the 
major  components  were  available  from  the  laboratory  but  the  Interfacing 
and  peripheral  components  required  to  complete  the  system  were  either 
designed  and  fabricated,  or  purchased.  Software  for  aPDP  11/40  digital  < 

computer  Interfaced  with  a fringe  reader  and  analyzer  was  developed  and 
utilized  to  digitize  the  film  and  fringe  patterns,  plot  the  fringe  density, 
locate  the  flaws  on  the  plot,  determine  the  flaw  size,  and  compute  the  j 

surface  strains. 

In  summary,  fringe  data  were  taken  using  both  double-exposure  and  \ 

real-time  optical  holography.  Young's  Fringe,  and  shearing-speckle  j 

Interferometry.  The  data  were  reduced  by  Computer  and  the  flaws  were  j 

detected,  located  and  the  approximate  flaw  size  determined.  Strains, 
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displacements,  or  stresses  were  determined.  Optical  holography  yielded 
the  best  fringe  resolution  and  the  speckle  techniques  simplified  the  data 
reduction  to  strain. 

For  application  at  the  end  of  the  production  line,  the  speckle 
methods  are  the  most  adaptable  because  they  can  be  more  easily  digitized 
using  a rotating  mirror  system  or  a vldicon  system  with  a scan  converter 
interfaced  with  a computer.  As  the  system  now  exists,  the  shearing 
speckle  techniques  are  used  to  Inspect  the  whole  surface  for  flaws. 

It  is  expected  that  907.  to  957.  of  the  items  will  not  contain  flaws  and 
can  be  accepted  with  no  further  testing  for  flaws.  The  remAining  57.  to 
107.  of  the  items  are  then  analyzed  only  in  the  regions  where  flaws  are 
detected  using  Young's  Fringe  methods  to  compute  the  surface  strains 
and  accept  or  reject  based  on  a criterion  established  by  testing  some 
of  the  first  flawed  cylinders  from  the  production  line.  Data  from  these 
tests  are  more  reliable  because  of  the  variations  in  composite  materials 
resulting  during  fabrication  on  different  production  lines.  A data  bank 
from  these  tests  is  then  established  in  the  computer  memory  and  a compari- 
son of  the  surface  strains  versus  acceptable  strains  is  made  by  the  com- 
puter on  all  future  items  and  the  item  will  be  shown  as  accepted  or 
rejected  on  a printout. 

A system  to  be  used  at  the  end  of  a production  line  would  be  better 
optimized,  less  expensive, and  simpler  than  some  of  the  laboratory  compo- 
nents used  here.  Complete  documentation  and  drawings  required  to  purchase 
and  fabricate  any  cosiponents  are  on  file  in  this  laboratory  and  can  be 
made  available  to  personnel  in  Government  or  Government  contractors  in 
industry  for  the  fabrication  of  systems  to  be  used  in  production  line 
applications.  Procedures  for  operation  of  the  system  are  presented  in 
this  report. 

X.  FUTURE  SYSTEM  IMPROVEMENTS 

Further  refinements  of  the  system  should  be  directed  toward  the 
film  recording  to  speed  up  the  process.  It  is  possible  to  eliminate  the 
Young's  Fringe  analysis  and  use  the  shearing  speckle  interference  fringes 
to  do  both  inspection  and  analysis. 

It  may  be  possible  to  eliminate  the  film  completely  and  speed  up 
the  process  considerably  by  forming  the  interference  fringes  electronically. 
This  was  explored  by  forming  the  fringes  in  the  scan  converter  without 
success  but  it  is  still  considered  feasible  if  the  patterns  can  be  mixed 
on  the  target  of  the  vidlcon-silicon  tube  of  the  camera. 

The  capabilities  of  this  system  can  be  directed  toward  a production 
line  Inspection  technique  to  assure  that  the  design  material  properties 
of  composite  structures  are  obtained  in  all  the  items.  With  further 
development  and  refinement  of  the  techniques,  it. is  considered  feasible 
to  coiid>ine  the  flaw  detection  and  materials  into  one  test. 
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Appendix  A.  COMPUTER  CODE  FOR  VIDICON  CAMERA 


The  vldlcon  computer  code  for  Young's  Fr'inge  analysis  of  interfero- 
grams  Is  used  to  control  a vldlcon  camera  and  scan  converter  for  digitizing 
and  processing  video  Images.  The  program  controls  the  positioning  of  an 
X-Y  translation  table  and  Initiates  the  storage  of  video  Images.  The 
Images  which  are  stored  In  a scan  converter  are  digitized  and  then  maxima- 
minima  data  are  extracted  to  generate  fringe  spacing  Information  and 
localized  displacement  data.  The  flow  chart  for  the  program  Is  shown  In 
Figure  A-1.  Operation  of  the  computer  code  Is  as  follows.  After  the 
necessary  parameters  concerning  the  scan  boundary,  axis  step  Increments, 
and  X-Y  stage  translation  rates  are  given,  the  computer  requests  the 
locations  of  two  vertical  and  two  horizontal  scan  lines  for  the  extrac- 
tion of  video  data  from  the  scan  converter.  Next  the  machine  requests 
data  concerning  the  erasure,  storage,  and  filtering  of  video  data.  The 
final  piece  of  Information  to  be  provided  Is  the  system  scale  factor. 

Once  this  Information  Is  given,  the  computer  computes  the  Initial  coordi- 
nates for  starting  the  raster  scan  of  the  Inter ferogram.  It  then  advances 
the  X-Y  stage  to  the  calculated  coordinates.  The  Interferogram  aperture 
is  4 X 5 In,  and  the  total  number  of  available  points  for  scanning  in 
the  X and  Y directions  Is  4000  and  5000,  respectively.  This  gives  20 
million  possible  data  points  for  scanning. 

After  positioning  the  X-Y  Table  for  each  scan,  the  computer  goes 
through  a sequence  of  operations  for  extracting  data.  First,  the  scan 
conversion  memory  Is  erased  and  then  one  frame  of  video  from  the  vldlcon 
camera  Is  stored.  Next,  the  computer  scans  along  two  horizontal  and  two 
vertical  lines  In  the  scan  converter  extracting  the  video  levels  stored 
there.  The  data  for  each  of  the  two  sets  of  scan  lines  Is  reduced  to  a 
zero  base  level,  averaged  and  squared.  The  data  signals  for  the  hori- 
zontal and  vertical  directions  are  then  filtered  and  the  locations  of 
the  fringe  maxima  and  minima  are  determined.  From  this  Information,  the 
fringe  spacing  and  orientation  are  computed  and  deflection  data  are 
generated.  The  results  are  printed  out  and  the  X-Y  stage  advances  to 
a new  location  for  scanning.  All  scanning  of  the  Interferograms  Is  per- 
formed In  raster  fashion  as  shown  In  Figure  A-2.  After  a raster  scan  Is 
completed,  the  computer  returns  to  Its  original  position.  An  Important 
point  to  note  Is  that  the  location  of  the  X-Y  stage  before  the  raster 
scanning  begins  Is  considered  to  be  the  origin.  This  Increases  the  flexi- 
bility In  using  the  scanner. 

Another  major  point  In  the  operation  of  extracting  fringe  data  Is 
the  filtering  technique  which  Is  used.  The  Input  signal  from  the  scan  i 

converter  to  the  ARlli^)  converter  Is  noisy  due  to  speckle  noise  and  some  elec-  i 

trlcal*  noise  . Therefore,  the  digitized  scan  data  must  be  filtered . The  i 

filtering  Is  accomplished  by  specifying  the  amplitude  of  the  noise  which 

Is  superimposed  upon  the  desired  signal.  This  parameter  Is  referred  to  j 

as  the  filter  level  and  may  be  determined  from  the  relation: 


Filter  level  ■ noise  level  (V)  * ^ — 
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Figure  A-1.  Vidicon  computer  code  flow  chart  for  Young's  Fringe 
analysis  of  interferograms. 


TT:=DX1 JVTEST.FOR 


C VIDICON  YOUNG'S  FRINGE  ANALYZER 

DIMENSION  ID<2fl024)rIDD(2rl024) 

WRITE(5>1) 

1 FORMATC'  GROUND  EQUIPMENT  AND  MATERIALS  DIRECTORATE') 

WRITE(S>2) 

2 FORHATK'  VIDICON  YOUNGS  FRINGE  ANALYZER  SYSTEM') 

URITE<5>3) 

3 FORMAT ('  ') 

WRlTE<5f4) 

4 FORMAT('  INPUT  THE  FOLLOWING  SCAN  PARAMETERS:') 

URITE(5t3) 

URITE(5»5) 

5 FORMAT('  X-AXIS  BOUNDARIESr  Y-AXIS  B0UNDARIES-4I5: ' ) 

READ(5.6)  IX1»IX2.IY1»IY2 

6 FORMAT < 415) 

URITE<5f7) 

7 FORMAT('  X*Y  AXIS  INCREMENTS*  X»Y  AXIS  STAGE  RATES-4I5:') 
READ(5*A)  IXCf lYCf IXR*IYR 

URITE<5f8) 

8 FORMAT<'  XY-1  AND  XY-2  SCAN  INTERSECTS-4I5; ' ) 

READ(5f6)  IXX1.IYY1.IXX2*IYY2 

WRITE<5r9) 

9 format:'  store  delay*  prime  delay*  no.  primes*  filter  LEVEL-415 
READ<5*6)  ISD*IPD*IPC*IT 

WRITE(5*18) 

18  format:'  system  scale  FACT0R-F10.4:') 

READ:5*19)  sf 

19  F0RMAT:F10.4) 

C MOVE  TO  THE  INITIAL  START  POSITION 

M1=IX1-IXC 
M2=IY1-IYC 
CALL  xadv:mi*ixr) 

CALL  YADV:M2*IYR) 

C PERFORM  THE  INTERFEROGRAM  SCAN 

M3=::iY2-IYl)/IYC)+l 

M4=l 

M5=::ix2-ixi)/ixc)+i 

M6=l 

DO  22  J°1*M3*1 
CALL  YADV:iYC*IYR) 

DO  21  I=M4*M5*M6 
CALL  XADV:iXC*IXR) 

26  CONTINUE 

CALL  store:isd*ipd*ipc) 

CALL  scan:ixxi*iyyi*id) 

CALL  SCAN:IXX2*IYY2*IDD) 

IXMIN-1024 

IYMIN=1024 

IXXMINS1024 

IYYMIN=1024 

DO  10  KK=1* 1024*1 

if:ixmin.gt.id:i*kk))  ixmin=id:i*kk) 
if:iymin.gt.id:2*kk))  iymin=id:2*kk) 
if:ixxmin.ot.idd:i*kk)>  ixxnin=idd:i*kk) 

10  if:iyymin.qt.idd:2*kk))  iyymin“Idd:2*kk) 

DO  11  KKK>1* 1024*1 

id:i*kkk)>*id:i*kkk)-ixmin 

id:2*kkk)-id:2*kkk)-iymin 

idd:i>kkk)’=idd:i*kkk)-ixxmin 

11  id:2*kkk>-'Idd:2*kkk>-iyymin 


DO  12  K>lfl024fl 
ID(lfK>*><ID(lfK)+IDD<l»K)  >/2 
ID(2fK)=<ID(2»K)+IDD(2»K))/2 
ID<lrK>>ID<lrK)»ID(lfK) 

12  ID(2fK>=ID<2fK)«ID(2FK> 

IXNAX=0 
lYNAX^O 
ICX=1 
ICY«1 

DO  13  KKKK<lrl024f 1 
IFdXHAX.LT.IDdrKKKK))  ICX^^KKKK 
IFdXMAX.LT.IDdfKKKK) ) IXMAX=IDd  fKKKK) 

IF < I YMAX . LT . ID < 2 F KKKK  > ) ICY=KKKK 

13  IFdYMAX.LT.ID(2FKKKK))  IYMAX=ID<2fKKKK) 

CALL  PROCESdDFlDOFlTFNNXFNNY) 

NDX>102A 

IF<NNX.LT.3)  GOTO  15 
NNXM=NNX-1 
DO  14  II-2fNNXMf1 

14  IFdDDdFiD.EQ.ICX)  NDX=IABSdDDd  f I I + l >-IDDd  f II-l  > ) 

15  NDY>1024 
IF<NNY.LT.3)  GOTO  17 
NNYM-NNY-1 
DO  16  JJff2fNNYMf1 

16  IF(I0D(2FjJ>.Ea.ICY)  NDY=IABS( IDD< 2f Il-d  )-IDDd f II-l ) > 

17  CONTINUE 
DX»SF/FLOAT(NDX) 

DY=SF/FLOAT<NDY) 

NX=IXl  + <I-l))|cIXC 
NY=>IY1  + (J-1)*IYC 
URITE(5f20)  IfJfDXfDYfNXfNY 

20  F0RMAT(3H  X=fI5f3H  Y=fI5f4H  DX=fE12.5f 

14H  DY>fE12.5f4H  XS>fI5f4H  YS=fI5) 

ITEST»IPEEK<  *177570) 

IF(ITEST.NE.l)  GOTO  25 
URITE(5f23) 

23  FORMAT('  NEW  FILTER  LEVEL:') 

READ(5f24)  it 

24  FORMAT (IS) 

GOTO  26 

25  IF(ITEST.EQ.2)  GOTO  26 

21  CONTINUE 
CALL  XADVCIXCfIXR) 

M7-M4 

M4*M5  J 

M5-M7  j 

M6-(-l)»M6  j 

22  IXC«(-1)*IXC  ! 

M2— NY  i 

Ml— NX+IXC  j 

CALL  XADV(MIfIXR) 

CALL  YADV(M2fIYR) 

STOP 

END 

« 


V 


fM  nn 


tt:=dxi:xadv.for 


SUBROUTINE  XADUdSrIR) 
IS=NO.  STEPS  (+=FWD»-=REV) 
IR-ADVANCE  RATE  OF  STAGE 
X=0. 

IFCIS.GT.O)  GOTO  3 
IP>‘IABS(IS> 

DO  2 I>lfIPfl 

CALL  IPOKE< •167772F '000000) 
CALL  IPOKE< ■167772f *020000) 
DO  1 J^lrlRrl 
Y=SIN<X) 

CONTINUE 
GOTO  6 ■ 

CONTINUE 

DO  S Il^lflSfl 

CALL  IPOKEC ■167772» '010000) 

CALL  IPOKE< '167772f *030000) 

DO  4 JJ=1pIRf1 

Y«SIN<X) 

CONTINUE 

CONTINUE 

RETURN 

END 
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tt:-dxi:yad‘'.for 


I 


f 


c- 

c 


1 

2 

3 


4 

5 

6 


* 


SUBROUTINE  YADUdStlR) 
■IS»NO.  STEPS  <+-FUD»-»REU> 
■IR>AOUANCE  RATE  OF  STAGE 
X-0. 

IFdS.GT.O)  GOTO  3 
IP-IABS(IS) 

DO  2 I-lfIPfl 

CALL  IPOKr(  *167772. ■000000). 
CALL  IPOKE( *167772. *100000) 
DO  1 J>l.IR.l 
Y=SIN<X) 

CONTINUE 

GOTO  6 

CONTINUE 

DO  5 11=1. IS. 1 

CALL  IPOKE( *167772. *040000) 

CALL  IPOKEC 167772. *140000) 

DO  4 JJ=l.IR.l 

Y=SIN(X) 

CONTINUE 

CONTINUE 

RETURN 

END 
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TT:=DX1 JSTORE.FOR 


SUBROUTINE  STORE ( ISO r IPO > IPC) 

DO  1 In^lrlPCfl 

CALL  IP0KEC167772»'00) 

CALL  IPOKEC 167772 t'OA) 

DO  1 J=l>IPDfl 
Y=0. 

1 X=SIN<Y) 

CALL  IPOKEC ■167772f *00) 

CALL  IP0KEC167772.*01> 

DO  2 K=lrISDfl 
Y=0. 

2 X=SIN<Y) 

CALL  IPOKEC ■167772F '00) 

RETURN 

END 

* 


[ 

I TTJ=DXi:SCAN.FOR 


SUBROUTINE  SCANC IX> lYr ID) 

DIMENSION  IDC2fl024) 

C X SCAN  THE  SCAN  CONVERTER 

CALL  IPOKEC •170414»IY) 

DO  1 I>FlFl024rl 
CALL  IPOKEC •1704I2F I) 

CALL  IPOKEC ■167772f ’10) 

CALL  IPOKEC* 147772* '30) 

CALL  IPOKEC *170400* *20401) 

IDC1*I)=IPEEKC *170402) 

1 CALL  IPOKEC *147772* *00) 

C Y SCAN  THE  SCAN  CONVERTER 

CALL  IPOKEC *170412*IX) 

DO  2 J*=l*  1024*1 

CALL  IPOKEC *170414* J)  j 

CALL  IP0KEC*147772**10) 

CALL  IPOKEC *147772* *30) 

CALL  IPOKEC *170400* *20401) 

IDC2*J)=IPEEKC *170402) 

2 CALL  IPOKEC *147772* *00) 

RETURN 
END 

« 
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tt:=dxi{proces.for 


SUBROUTINE  PROCES ( ID f I DD 1 1 T f NNX r NNY ) 
DIMENSION  ID(2f 1024) rIDD<2f 1024) 

K=I 

1 III=0 
LC»I 
LCl  = l 
IDI=0 
J=1 

2 J=J+1 

ID2=IABS<ID(K> J>-ID(KfLC> ) 

IF<ID2.GT.ID1>  LC1=J 
IF<ID2.GT.ID1>  ID1=ID2 
ITP=IABS(ID<K.J)-ID<K»LC1)) 

IF(ITP.LT.IT)  GOTO  3 
III=III+1 

IF ( ID(K>LC1 > .LT.ID<K»LC) ) IDD<K» III >=-LCl 

IF(ID<KfLCl > .GT.ID(KfLC) ) IDD(K r I I I )=LC1 

LC=LC1 

ID1=0 

J=LC 

3 IF( J.LT.1024)  GOTO  2 
IF<K.EQ.l)  NNX=III 
NNY=III 

K=K+1 

IF(K.LE.2)  GOTO  1 

RETURN 

END 
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Appendix  B.  COMPUTER  CODE  FOR  ROTATING  MIRROR 


The  computer  code  for  analyzing  data  from  the  rotating  mirror  assem- 
bly converts  speckle  fringe  data  into  displacement  Information.  The 
operation  of  the  code  is  described  as  follows.  Before  an  interferogram 
scan  can  begin,  the  operator  of  the  system  is  requested  to  input  data 
parameters.  The  computer  code  scans  in  a raster  type  of  operation  in 
which  up  to  20  million  data  points  may  be  examined.  Initially,  the 
operator  of  the  computer  program  sets  up  the  boundaries  for  the  scan  by 
giving  the  computer  the  left  most,  right  most,  lower,  and  upper  coordi- 
nates. Next,  the  Increments  between  scan  points  in  the  X and  Y direc- 
tions are  given  along  with  the  rates  for  translating  the  X and  Y stages. 
This  information  instructs  the  nuchine  for  the  spacing  between  points 
for  which  data  are  desired  and  how  fast  the  operation  of  translating  the 
X-Y  stages  is  to  occur.  Finally  the  computer  requests  a scan  level  and 
system  scale  factor.  Once  the  necessary  information  has  been  provided, 
the  computer  calculates  the  initial  coordinates  for  scanning  and  advances 
to  them. 

A sequence  of  operations  is  required  to  raster  scan  an  interferogram. 
The  flowchart  for  the  program  is  shown  in  Figure  B-1  and  Figure  B-2.  The 
Important  point  to  note  from  the  figures  is  the  method  in  which  fringe 
maxima  are  located.  Once  the  X-Y  Table  has  been  properly  positioned, 
the  a/D  converter  starts  to  sample  a waveform  similar  to  the  one  shown 
in  Figure  B-3.  First  the  program  for  converting  the  waveform  remains 
inactive  until  the  signal  voltage  swings  from  below  the  scan  level  I to 

above  it.  The  A/D  converter  then  samples  the  waveform  for  1000  conversions . 
Finally  the  program  converts  the  sampled  data  into  fringe  spacing  by  locating  the 
maxima  and  computing  the  spacing.  The  criteria  for  determining  these 
qiaxlma  are  indicated  in  Figure  B-4.  Whenever  the  slope  between  three 
successive  data  point  samples  changes  from  positive  to  negative,  a maxima 
exists.  This  technique  for  maxima  determination  requires  that  the  input 
signal  be  filtered. 
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Figure  B-1.  Rotating  mirror  analyzer  computer  code  flow  chart. 
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rNOTOOETECTOR  VOLTAGE  V 


C ROTATING  HIRROR  YOUNG'S  FRINGE  ANALYZER 

URITE<Srl) 

1 FORHAT('  GROUND  EQUIPMENT  AND  MATERIALS  DIRECTORATE') 
URITE<Sr2> 

2 FORMAT<'  ROTATING  MIRROR  YOUNGS  FRINGE  ANALYZER  SYSTEM') 
URITE(Sf3) 

3 FORMAT ('  ') 

URITE(S>4) 

4 FORMAT<'  INPUT  THE  FOLLOWING  SCAN  PARAMETERS:') 

WRITE(5>3) 

WRITE(SfS) 

5 format:'  X-AXIS  BOUNDARIESf  Y-AXIS  B0UNDARIES-4I5 I ' ) 
READ<S>A)  IXl>IX2rIYlf IY2 

6 format: 415) 

WRITE:Sf7) 

7 format:'  X*Y  axis  increments*  X*Y  axis  stage  RATES-415:') 
READ:Sr6)  IXC> lYC* IXR* lYR 

urite:s>8) 

8 format:'  scan  level*  system  scale  FACT0R-IS*F10. ' ) 

REA0:5*9)  IS*SF 

9 F0RMAT:IS*F10.4) 

C MOVE  TO  THE  INITIAL  START  POSITION 

Ml-IXl-IXC 
M2-IY1-IYC 
CALL  xadv:mi*ixr) 

CALL  YADV:M2*IYR) 

C PERFORM  THE  INTERFEROGRAM  SCAN 

M3*: :iY2-IYl)/IYC)+l 

M4-1 

M5*::ix2-ixi)/ixc)+x 

MA-1 

DO  22  J-1*M3*1 
CALL  YADV:iYC*IYR) 

DO  21  I*M4*MS*M6 
CALL  XAOV:iXC*IXR) 

26  CONTINUE 

CALL  rload:d*is) 

D-SF/D 

NX-IXlf :i-l)»IXC 

ny-iyi+:j-1)*iyc 
write:s*20)  i*j*nx*ny*d 

20  F0RMAT:3H  X-*I5*3H  Y**IS*4H  XS**IS* 

14H  YS-*IS*3H  D**E12.S) 

ITEST-IPEEK: *177570) 

IF:ITE8T.NE.2)  goto  25 
HRITE:5*23> 

23  format:'  new  scan  LEVEL}') 

, READ:5*24>  is 

I 24  format: 15) 

(GOTO  26 

25  IF:ITE8T.EQ.1)  goto  26 

21  CONTINUE 

CALL  XADV:iXC*IXR) 
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I 


H7-H4 

M4-H5 

M5-H7 

22  IXC>>(-1>«IXC 

M2—MY 
HI— NX+IXC 
CALL  XADV(HltlXR) 
CALL  YADV(H2tIYR) 
STOP 
END 

« 


tt:=dxi:rload.for 


SUBROUTINE  RLOAD(DflS> 

DIMENSION  ID(1000>rIL(32) 

7 CALL  IPOKE< •170400f *21001) 

IT=IPEEK( *170402) 

IF(IT.GT.IS)  GOTO  7 

2 CALL  IP0KE(*170400f *21001) 

IT»IPEEK< *170402) 

IFdT.LT.IS)  GOTO  2 

DO  3 I=l>1000rl 

CALL  IPOKE( *170400» *21001) 

3 ID(I)=IPEEK( *170402) 

DO  4 J=1>1000f1 

4 IF(ID<J).LT.IS)  ID(J)=0 
KC=0 

NC=0 

KT=0 

DO  5 K=2p1000f1 

IF(NC.GT.32)  GOTO  5 

IF(ID(K).GT.ID<K-1))  KC=1 

IF(ID<K) .LT.ID<K-1) .AND.KC.EQ.l)  KT=1 

IF(KT.NE.l)  GOTO  5 

NC=NC+1 

KC=0 

KT=0 

IL(NC)=K-1 

5 CONTINUE 
IF<NC.LT,2)  D=1E6 
IF(NC.LT,2)  GOTO  A 

D=FLOAT  <IL(NC)-IL(1)) /FLOAT ( NC-1 ) 

6 RETURN 
END 


on 


tt:=dxi:xadv.for 


SUBROUTINE  XADUdSfIR) 
-IS=NO.  STEPS  <+=FWD»-=REV> 
-IR= ADVANCE  RATE  OF  STAGE 
X=0. 

IF(IS.GT.O)  GOTO  3 
IP=IABS<IS) 

DO  2 I=1»IP»1 

CALL  IP0KEC1A7772»  *000000) 
CALL  IPOKE(  •l<S7772f  *020000) 
DO  1 J^lrlRd 
Y=SIN(X) 

CONTINUE 

GOTO  6 

CONTINUE 

DO  5 II=lfISfl 

CALL  IPOKE(  *l<S7772f  *010000) 

CALL  IPOKE( *167772» *030000) 

DO  4 JJ=ldR>l 

Y=SIN(X) 

CONTINUE 

CONTINUE 


1 


I U 


no 


1 


tt:=dxi:yadv.for 


SUBROUTINE  YADU<IS>IR> 

IS=NO.  STEPS  (+=FWDf-=REV) 

IR=ADVANCE  RATE  OF  STAGE 

X=0. 

IFdS.GT.O)  GOTO  3 
IP=IABS(IS) 

DO  2 I=lrIPrl 

CALL  IPOKE( •167772f *000000) 
CALL  IPOKE( •167772» *100000) 
DO  1 J=1«IR<1 
Y=SIN<X) 

CONTINUE 

GOTO  6 

CONTINUE 

DO  5 II=lfIS»l 

CALL  IPOKE( *167772» *040000) 

CALL  IPOKE( *167772f *140000) 

DO  4 JJ=l»IRfl 

Y=SIN(X) 

CONTINUE 
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4308  Governors  Drive 

Attn:  T.  Ward  1 

Huntsville,  Alabama  35805 

DRSMI-LP,  Mr.  Voigt  I 

DRDMI  -X  1 

-T,  Dr.  Kobler  1 

-TL,  Mr.  Lewis  1 

-TLA,  Mr.  Pettey  1 

-TLA,  Dr.  Mulllnlx  175 

Mr.  Schaeffel  12 
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